{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Ill-conditioning of Vandermonde matrix\n", "\n", "\n", "**Randall Romero Aguilar, PhD**\n", "\n", "This demo is based on the original Matlab demo accompanying the Computational Economics and Finance 2001 textbook by Mario Miranda and Paul Fackler.\n", "\n", "Original (Matlab) CompEcon file: **demlin02.m**\n", "\n", "Running this file requires the Python version of CompEcon. This can be installed with pip by running\n", "\n", " !pip install compecon --upgrade\n", "\n", "Last updated: 2022-Ago-19\n", "
\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "from numpy.linalg import norm, cond, solve\n", "import time\n", "import matplotlib.pyplot as plt\n", "np.set_printoptions(precision=4)\n", "plt.style.use('seaborn-dark')\n", "plt.style.use('seaborn-talk')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Compute approximation error and matrix condition number" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "nvalues = np.arange(6, 101)\n", "\n", "data = pd.DataFrame(index=nvalues, columns=['conv', 'errv'], dtype=float)\n", "data.index.name = 'n'\n", "\n", "for n in nvalues:\n", " v = np.vander(1.0 + np.arange(n))\n", " data.loc[n, 'errv'] = np.log10(norm(np.identity(n) - solve(v, v), ord=np.inf))\n", " data.loc[n, 'conv'] = np.log10(cond(v, p=2))\n", "\n", "data.rename(columns=dict(errv=\"$||I-V^{-1}V||_\\infty$\", conv='Condition number of $V$'), inplace=True) " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Smooth using quadratic function" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Users\\randa\\AppData\\Local\\Temp\\ipykernel_18304\\1182016897.py:2: FutureWarning: `rcond` parameter will change to the default of machine precision times ``max(M, N)`` where M and N are the input matrix dimensions.\n", "To use the future default and silence this warning we advise to pass `rcond=None`, to keep using the old, explicitly pass `rcond=-1`.\n", " b = np.linalg.lstsq(X, data)[0]\n" ] } ], "source": [ "X = np.vstack([np.ones_like(nvalues), nvalues, nvalues**2]).T\n", "b = np.linalg.lstsq(X, data)[0]\n", "data_smooth = pd.DataFrame(10 ** (X @ b), index=data.index, columns=data.columns)\n", "\n", "#b = np.linalg.lstsq(X, conv)[0]\n", "#conv = 10 ** (X @ b)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot matrix condition numbers" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABAYAAAHnCAYAAAArel9ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/NK7nSAAAACXBIWXMAAA9hAAAPYQGoP6dpAACmqElEQVR4nOzdd3jb5bk+8Ftb3ns73o5jO45tyCBkL2c6TgizUMooqwdooZRDKWV0wKEt0DLKrwFKWmgZSSDbSZxN9p7ee+89tL+/P2SLBEmO7djy0P25rnOd4u9j6bEsO9at931ekSAIAoiIiIiIiIjILomHuwEiIiIiIiIiGj4MBoiIiIiIiIjsGIMBIiIiIiIiIjvGYICIiIiIiIjIjjEYICIiIiIiIrJjDAaIiIiIiIiI7BiDASIiIiIiIiI7xmCAiIiIiIiIyI4xGCAiIiIiIiKyYwwGiIiIiIiIiOwYgwEiIiIiIiIiO8ZggIiIxhxBEIa7BSIa4fh7gojoewwGiIhGoJiYGMTExFzzsRdeeAExMTHYuXOnxf/uj66uLmzcuBEPPPAA5s6di4kTJ2L69Ol49NFHsXv37kH5GgbTJ598gpiYGLz33numj1n7+jdv3oyXXnrpmo/dyGM1UMNxn4MpLy8P999/P5KTk5GcnIwXXnjBYt3DDz+MmJgYfPjhh9e9za6uLiQnJyMuLg41NTWD3XK/WHpOkbmhfJy++eYb0++6JUuWXLf+rbfeMtV/8803N3Tfln5P9Gaofp7Pnj2LmJgYLFq0qE/1f/nLXxATE4M77riDwQYRDSoGA0REdiYzMxMrV67Eiy++iEuXLiEwMBDz5s1DYGAgDh06hKeeegrPPvss9Hr9cLfab6dOncLzzz+P+vr64W5l1Hv66adx4sQJBAUFYfbs2UhISLBYd9tttwEAtm/fft3bzMjIQGdnJ2bOnAk/P79B7ZdGt6KiIuTl5fVaM1gvykfS74nx48cDAMrLy6FSqXqtraysxL///W8AwIsvvgiRSDTk/RGR/ZAOdwNERGQ7BQUF+NGPfoSuri489thjePTRR+Hs7Gy6np2djaeeegrbt2+Hs7Mzfve73w1jt7179tln8cgjj1zzAtNgMPS5dqgNx30Olvb2dhQWFsLR0REbN26EQqGwWrtw4UK4uroiLy8POTk5ZitdrrZlyxYA34cJRADg6uqK1tZW7Nq1C9HR0RZrrly5gtLSUshkMmi12hu6P2u/J3ozVD/Pzs7OCAwMRGVlJQoKChAfH2+19q233oJarcaKFSuQnJw8qH0QEXHFABGRnRAEAc899xy6urrw9NNP49lnn70mFACACRMm4KOPPoJMJsP69etRVFQ0TN1en6+vLyIjI82+hhutHSzDcZ+DRaPRAAA8PDx6DQUAQKFQYOnSpQB6XzXQ0NCAY8eOwd3dHfPnzx+8ZmnUmz17NqRSaa/bmHpWC8yYMcNWbV1jKH+ee8KQ3lZMXLp0Cdu3b4dSqcRzzz036D0QETEYICKyE6dOnUJmZib8/f3xyCOPWK0LCwvD3XffjTlz5qCuru6aawaDAV999RXuuOMOJCUlITk5GXfeeSfWr19vtt+1Z09ucXEx/vOf/yA1NRUJCQmYOXMmXn75ZTQ0NJjdd1dXF9577z2kpKRg0qRJSE1NxebNmy32aWnmwv333w8AOHDgAGJiYkz74q3tDx7Kr8fSfQ7kMdFqtfjkk0+wfPlyJCYmYv78+Xj33XdRWlp6zdd4PX39Wn/84x9j+vTpAICKigqL8y5+aM2aNQB6Dwa2bdsGnU6H1NRUyOVy08cbGxvx1ltvITU1FcnJyUhISMD8+fPx29/+FlVVVdfcxlA/pwBjgLZx40bcddddpvkK99xzD3bs2GFW29PPuXPn8PDDD5t62bdv3zXXy8rKsG7dOixZsgQJCQlYsGAB1q5dC0EQ0NzcjJdeegnTp0/H5MmTcf/99+PKlStm9zWSfvb6+zhdj5ubG2655Rbk5OSguLjYYs3OnTsxadIkBAUFWbze1+dRX35PWPp+/vDn+cKFC4iLi8PEiRPNXtCfPn0asbGxmDp1ap9mafRsJygoKLBa8+abb0IQBDz00EMICAi47m0SEfUXtxIQEdmJ9PR0AMal31e/MLPE0lAuvV6Pp59+Gnv27IGjoyOmTZsGADh58iReeuklHDlyBG+//TbE4msz5zfffBP79u3DpEmTMGfOHJw8eRJfffUVzp07h2+//RZSqfGfIrVajYceeghnz56Fj48P5s6di7KyMjz//POIioq67teXnJyMmpoaHD16FH5+fpg6dWqvy22H+uvpTV9vQ6/X46mnnsL+/fvh7u6O2bNno66uDh988AG+++67697PQL7WW2+9FR4eHti1axccHR2xYMGC695+YmIiIiMjUVBQgAsXLiAxMdGspudFZk+IAAA1NTW46667UFVVhcjISMycORPt7e24cOECvv76axw6dMi0rWUgj19/n1M9q2q2bdsGJycnJCcnQyqV4uTJk3jmmWdw+fJlPP/882af98ILL6CtrQ2zZ89GZmam2XLw3/3udzh8+DCmTJmCoKAgHD9+HG+99RZUKhXS09PR1NSEpKQklJeX48SJE7jvvvuQnp4Of3//fn//RvLj1JvFixfj8OHD2L17Nx599NFrrvVsI7jnnntQXl5u9rn9eR715feEpe/nD1czJCYm4ic/+Qn++c9/4rXXXsPnn38OwBiwvPjiizAYDHjllVf6tPXgeisGMjIycOrUKfj6+vYa6hIR3QgGA0REdqJnW4C1IXLX869//Qt79uxBTEwMPv74Y/j6+gIw/lH+8MMPIz09HcnJyfjJT35yzecdPnwYH3/8MWbNmgUAqK+vx2233Ybc3FwcOXIEc+bMAQCsW7cOZ8+exa233oq///3vcHBwAGCciv6nP/3puv3dddddCAsLw9GjRxEbG4u//OUvw/r19Kavt7F+/Xrs378fSUlJ+Oijj+Dq6goA2LVrF5555pnr3s9AvtYnnngCjY2N2LVrFzw8PK77OPZYtWoV3nrrLWzbts0sGCgsLMSVK1cQGxuL2NhY08fff/99VFVV4cknn8RTTz1l+nhjYyPuvfdeFBYWYt++fVi5cuWAHr/+Pqe++OILbNu2DUlJSXj//ffh4+MDAKiursZDDz2ETz75BNOmTTP7Hre0tGDbtm3w9vaGwWAwe4F+7NgxfPbZZ5g8eTIAYOPGjXjxxRfxwQcfYNKkSfjyyy/h5uYGg8GAxx57zPRC9uGHHwYw8n72Bvo49WbhwoV49dVXLQYD6enpEIlEWLJkCT7++GOzz+3P86gvvyeu9/3s8fOf/xz79u3DqVOnsGnTJtPPQElJCZYuXYrly5f36WvvWZGTn59vdk2r1Zp6/OUvfwlHR8c+3SYRUX9xKwERkZ2ora0FAHh5eQ3o8z/77DMAwJ/+9CfTCxMA8PPzw5tvvgkA+PTTT80+b+XKlaYXJgDg7e2NxYsXA7j2D+Gvv/4aIpEIv//9700vTADjcXhDMWhrqL+e3vT1NnrehXzjjTdMoQBgfHf16nfer2egX2t/rFq1ChKJBOnp6WbD3XpWC/xw6KCnpyfmzJlj9i6op6en6fH44XYCYOieU+vWrYNIJMKf/vQn04tdAPD39zetovnXv/5l9nmLFy+Gt7c3AFh8EZmammoKBQBg2bJlpv/97LPPws3NzfS5CxcuBACUlJSYakbaz95AH6feeHp6YsqUKbh06RIqKyuvubZr1y4kJiYiMDDQ6ucO5HlkzfW+nz2USiX++Mc/QiQS4c9//jMOHjyI//znP/Dx8cGrr77a5/uLiIiAVCpFeXk5urq6rrn2xRdfoLi4GBMnTkRaWlqfb5OIqL8YDBAR2YmeZcM6na7fn1tZWYnKykqEhoZiwoQJZtfj4+MRGhqKqqoqs6W+lpaV97yY6PkjuLq6GuXl5YiMjERwcLBZfV+Ws/fHUH8919OX26ivr0deXh4iIiIQERFhVt/zgud6buRr7Q9fX1/MmDEDdXV1OHHihOnjgiBg27ZtkMlkSE1NveZznnnmGaxduxZKpdL0scbGRnz33XemffaWJtAPxXOqpqYGJSUlCAgIQGhoqFn9lClTIJfLcebMGbOjPK9eBWHJD1fpODg4mF6A//B70hMA9QyAHGk/ezfyOF3PkiVLAOCaZfuXL19GaWmpacClJQN9Hllzve/n1SZPnoz77rsP9fX1eOKJJ2AwGPCHP/wB7u7ufb4NuVyO0NBQCIKAwsJC08dbW1vxwQcfAODxhEQ09LiVgIhojDh9+jS+/PJLs4+npKQgJSUFPj4+yMnJQWNjY79vu2cIobXBXz3XSkpKUF9ff80LDBcXF7NaiUQCAKahaT2rGaztx7X0guVGDPXXcz19uY3q6moAsPouqbWP/9CNfK39tXr1ahw6dAg7duwwDTA8c+YMysvLsXjxYnh4eJh9TnFxMT7//HOcP38eJSUlaG1tBQDTiyBLj+lQPKd6hsRVVlZed9hiS0sLPD09Tf999WoOS3pWBFyt5+v74WPywxd/I+1n70Yep+tJSUnB7373O+zevRsPPPAAAOPQwZ5tBL0ZyPPImut9P3/o2Wefxc6dO1FXV4eFCxdi7ty5/fp8wDhnoKCgAHl5eaYZFR9++CGam5uxbNky3Hzzzf2+TSKi/mAwQEQ0RpSWlmLr1q1mHw8NDUVKSgri4uJw+PBhXLx48brL0C9duoSjR49i1qxZiIuLM/1R3ds7Vj01PxxsOBjvcvVloF9/DPfX05fb6FnZYe3M9b6+0LmRr7W/Fi5cCDc3N+zevRsvv/wyZDIZtmzZAgAWn3ObN2/Gr3/9a+j1ekRHR2POnDmIjo7GpEmTcP78efz1r3+1eD9D8ZzqeXfbx8cHt9xyS79uq7fl5sD3L8YHYrifq4P5OF2Pl5cXJk+ejNOnT6Ourg4+Pj7YtWsXkpKSTIMYLRno88ia630/fyg7O9t00sPx48dRXV3da7+WjB8/Hjt37jRt8SgrK8Pnn38OhULB4wmJyCYYDBARjRG33Xab2R7uq82bNw9r167FwYMHodFoen0R+MUXX2Djxo3Iz8/Hn//8Z9O+5t6WmvdcG8gMg54/oisqKixe73lXc7AM9dczGHoeE2t7o/tyDBpg269VLpdj2bJl+OKLL3DkyBHMmDEDO3fuhK+vL2bOnHlNbUdHB1577TXIZDJ8+umnpkn7PY4fP35DvfT3OdWzxN7d3b3PAxdtYaT97A3147R48WKcPHkSGRkZSExMRGlpKe677z6r9UP9PLoetVptOoVg8eLF2LVrF15++WWsXbu2X7fTc2RhTzDw1ltvQaPR4PHHH+91tQgR0WDhjAEiIjuRnJyMpKQkVFVV4Z///KfVuqysLNN59HfffTcA47L1wMBAlJSUIDs72+LnlJWVISQkpE/Hc/2Qr68vIiIiUFxcbHGA36FDh/p0O319h3Sov57B4O/vj5CQEBQVFVk8233//v19uh1bf6094VRGRgaOHTuGlpYW02DCqxUUFKCjowOJiYlmL+YEQcDRo0cBWF8xcT39fU4FBwfD398fxcXFKCsrM6svLi7GkiVL8NRTT/VrWfqNGmk/e0P9OKWkpEAsFmPPnj3YtWvXdbcRDOR5NJh79f/2t7+hqKgIK1aswDvvvIPY2FgcPHgQmzZt6tftXB0MnD9/Hunp6fDx8TE7oYGIaKgwGCAishMikci0vPudd97BX//6V7S3t19Tc+HCBTzxxBNQqVRYtWrVNftae961e+GFF1BfX2/6eG1tLf73f/8XAHDPPfcMuL/7778fAPDrX/8aLS0tpo9/8803OHDgQJ9uQ6FQAADa2tquWzvUX89g6OnxN7/5zTXfq6NHj+KLL74A0LcXObb8WidNmoTo6Gjs378fO3fuBGCcPfBDPS9is7KyTPMUAOPQvddffx0XL14EYHxHdqD6+5y67777oNVq8atf/eqaFRltbW349a9/jaKiIvj6+tp8CNxI+9kbysfJ19cXycnJOHnyJLZt24abbrqp18BjIM+j/vye6M3Fixexbt06uLq64te//jUkEgleffVViEQivP7666b5EH0xbtw4ODg4oLy8HL///e8BGGcXODk53VCPRER9xa0ERER2JD4+Hh999BH+53/+Bx9++CE+++wzTJw4Ee7u7igpKUFWVhYAYOnSpaY/Tns88MADOHPmDPbu3YtFixaZ3p07ceIEOjs7sWTJEtPAsIG46667cPjwYezZswcpKSmYOnUqamtrcf78eSQmJuLChQvXvY3g4GBIJBKcOXMGP/3pTzF58mQ8/vjjFmuH+usZDPfeey/27duH48ePY+HChZgyZQqamppw+vRpjBs3DqWlpX2av2Drr3XVqlX485//jM2bNyM5OdniqQp+fn6mpdfLli3D1KlTAQDnz59HU1MToqKikJ+ff80L4f7q73PqoYcewpkzZ7B//34sWbIEkyZNgoODA86cOYPW1lbEx8fjmWeeGXA/AzXSfvaG+nFavHgxzpw5g4qKCjz44IO91g7kedSf3xPWaDQavPjii9Dr9fjlL39pOt4wKSkJd955J7766iu8+uqrplMFrkcsFiMyMhKXL1/G5cuXER8fbzFQIyIaKlwxQERkZ6ZPn44dO3bgkUceQUhICC5evIiMjAzU1tZi3rx5+PDDD/HXv/7VbAaBRCLBe++9h1dffRURERE4fvw4Tp8+jdjYWLz55pv429/+1u+hXVcTi8V499138etf/xre3t44ePAgGhsb8Zvf/AYPPfRQn27Dy8sLr732GgIDA3H8+HHTMmJLhvrrGQxSqRRr167Fk08+CRcXF+zfvx+VlZV4+umnTQPJnJ2dr3s7tv5a09LSIJVKodPpeh10+X//93944okn4OPjg6NHjyIzMxNRUVF444038NVXX0EqleLQoUMDOmIT6P9zSiKR4IMPPsBrr72G6OhoXLx4ESdOnEBgYCB++ctf4vPPP+/T4z3YRtrP3lA/TosXL4ZIJIJYLO7TsZz9fR715/eENR988AHy8vKQnJyMu+6665prv/zlL+Hl5YU9e/Zgx44dfb7Nnu0EAI8nJCLbEwm23ChHREREfZadnQ1PT0/TALqrrVu3Dm+88QZeeeUV/OhHPxqG7oiIiGis4IoBIiKiEeq3v/0tZs2ahVOnTl3z8fLycqxbtw4ymQyzZs0apu6IiIhorOCKASIiohFqx44dePbZZyESiZCUlARfX180NTXh7Nmz0Ol0ePHFF02D44iIiIgGisEAERHRCHb69Gn8+9//xqVLl1BXVwc3NzckJCTgJz/5CaZPnz7c7REREdEYwGCAiIiIiIiIyI5xxgARERERERGRHWMwQERERERERGTHGAwQERERERER2THpcDdgL+rq2oa7BSIiIiIiIrIjPj4ufarjigEiIiIiIiIiO8ZggIiIiIiIiMiOMRggIiIiIiIismMMBoiIiIiIiIjsGIMBIiIiIiIiIjvGYICIiIiIiIjIjjEYICIiIiIiIrJjDAaIiIiIiIiI7BiDASIiIiIiIiI7xmCAiIiIiIiIyI4xGCAiIiIiIiKyYwwGiIiIiIiIiOwYgwEiIiIiIiIiO8ZggIiIiIiIiMiOMRggIiIiIiIi6gO9QUBeXTuau7TD3cqgkg53A0REREREREQjmUqrx9YrNfjiTDnKmlVwVUqx/sHJ8HSUD3drg4LBAA27zs5ObN36LTIydqG8vBQ6nQ5hYRFYsWIlVq68DWKx7Re23H57Kvz9A/D++2sBAE8++Siqq6uwYcNWU01TUyOUSgc4ODhYrbEnP3zMRiJBEPDhh+9h27bN0GjU+NnPfo7bbrtjuNsiIiIiohGqvkOD9ecqsPFCFVpUOtPHW1U6FNZ3wjOEwQDRDSstLcELLzyLysoKpKQsxdKlK6DVanD48CH85S//h3PnzuCVV/44LOHA1X7yk4fQ1aUy/fexY0fw2msv4dNP/2MKBn5YQyPP0aOH8d///hu33joTs2bNRWJiksW6jz/+f1i37mN88MFHSExMtlhjMBjw8MP3oaSkGJ999jWCgoKHsHMiIiIisqWC+g7890w50rNqodULZteXx/kiMch1GDobGgwGaNio1Wq8+OJzaG5uxscff4aoqGjTtbvvvg/vvfc2vvrqv4iJicOPfvTjYewUmDLllmv+OzPzMtrb23qtoZGnoCAPAPDYY08iMjLKal3PtcLCAqvBwPbtW5CXl4v77nuAoQARERHRGCAIAk6WNuM/p8txrLjJ7LpcIsLSOD/86OYgRHg5DUOHQ4fBAA2bb79dj+LiIvzmN69eEwr0ePzxp7B7905s2fIt7rnnPohEomHoksYSrdY4JMbR0bHXushI4/OxqKjA4vXOzg589NGH8Pb2wf33PzS4TRIRERGRTWn1BuzKrsV/z1Qgr67D7Lq7gwy3Jwbg9qRAeDmNja0DP8RggIbN3r274eDgiEWLlli8LpPJ8OGHn8DX1++aUODChXP45z8/QmbmZQACYmPj8dBDjyIp6SZTzZNPPgq5XIE777wHH330IYqKCuDu7oHly1fiwQcfuWZrwt69u/HZZ+tQWlqCoKAgPPPM82a9XD0/4I9/fBXp6dsAAHfcsRJJSTfh/ffXWpwxMNi9WuqrL59rbf+/pVkKjo6OSE1dhY8//gfKykoQGBiM//mfnyMxMRl///u72LcvAwAwbdp0PPPMr+Dq6nbNbW7btgn/+tenaGioQ2RkNH7608cxbdr0a2ouXjyPTz5Z2/24ABMnJuCRR55AXNzEa3qbMmUa9Ho99uzZDTc3N3zyyWfw9PSy+Fhc77G+/fZUVFdXmb5v/v4BVudBBAePg0KhQGGh5WBg3bpP0NjYgJdf/v11QwYiIiIiGplaurT45mIVvj5XifoOjdn1EA8H3HtzEJbF+UEpkwxDh7ZjF8cVHjt2DKtXr8ZNN92E5cuXY8+ePaZrFy9exOrVq5GUlIS77roLhYWFfbpmax0aHZo7tSPm/zo0uus33QtBEJCbm4OYmAmQSq3nU8HB4yCXf5/KHT58EE899Rhqaqpw//0P4YEHfora2hr8/OdP4PDhg9d8bmFhPl5++ddITr4ZP//5cwgMDMKnn36EjRu/NtXs2LEVr7zyIhQKBX72s6dw881T8fzzv0BjY6PVntLSbsPs2fMAAE8//Sx+8hPL7xgPdq/W3MjnWpKTk43XX/8d5syZh8cffxLNzU14+eUX8NxzT6OsrBSPPPIE5s1bgIyMnfjgg7/94HOz8M47f8aCBYvw6KM/Q3t7G55//hc4deqEqebEiWN4+unH0dnZjkceeRw/+clDqKmpxpNPPooLF85dc3t79uxCfn4enn76WaSmrrIaCvTlsX766V9e8317+ulfWn0MxGIxwsMjUVxs/jNfUVGODRu+REJCIlJSlvbtQSUiIiKiEaO8uQt/3puPFWtP4O+Hi81CgeRgN7y1Kh7rH5yM2xIDx3woANjBioGGhgY8/fTTePPNNzF37lwcO3YMTz75JDZv3gw/Pz/87Gc/w3PPPYdly5Zh7dq1eP7557Fhwwao1Wqr12zt//bk4ZsLVTAfeTF8RABuSwzACwvNtwD0RXNzM/R6Pby8vPv8OTqdDm+//Sd4e/vgk08+g5OTMwAgLW0N7r//Lrz11pu45ZYZpqChvr4O//d/b2PmzNkAgCVLlmPVqqXIyNiJO+64G3q9Hh9++B5iY+Pw/vtrIZPJAACxsXH4/e9fttrHxImTEBkZhUOH9mPWrLkICAgc8l57cyOfa0lDQz3efPMdzJgxCwAgkUjxzjt/glarxQcffGRavVFYmI8TJ45d87ldXV3485//iunTZwIAli1Lxd1334YPP3wPU6ZMg8FgwF/+8n+IjY3H+++vhURi/CW7Zs1dePDBH+Gvf/0zPv30v6bbU6vV+NOf3oG3t4/Vfvv6WM+ePRd5eTm9ft+uFhERiezsTDQ2NlwTSHzwwd+g0+nwi1/8qq8PKRERERENM0EQcLGyFZ+fLsfB/Aaz11YSEbAwxgc/ujkYcf4uw9LjcBrzKwYqKiqwdOlSzJ8/H2KxGDNmzEB4eDguX76M48ePw9nZGatWrYJcLscTTzyBkpIS5OTk9HrNlgRBwMYRFgoAgAAY+xIG1plEYnzq6fX6Pn9Obm42amtrsGbNnaYXfwDg7OyM2267A3V1tcjOzjR9XKlU4tZbZ5r+W6FQIDQ0DM3NTabba2pqxLJlK02hAAAsWrQE7u7uA/q6hqrX3tzI51oilyuuWfofEhIKAJg9e941WzoCA4PQ0FB/zeeGh0eYQgEAcHV1Q0rKUuTmZqOhoR65uTmoqqrArFlz0dbWhubmZjQ3N0OtVmPGjFnIy8tFbW2N6fODgoJ7DQWA/j/WfdUz9+Lq7QRnz57GoUP7sWJFGmJiJvT7NomIiIjItnQGARk5dXjoi/P46ZcXcOAHoYCTXIL7Jgdj00+n4g/LY+0yFADsYMXApEmTMGnSJNN/l5WVIT8/H9HR0Thy5AgiIiJM1yQSCcaNG4eCggLU1tZavRYTE2Oz/kUiEdYkBozIFQNrEgMGPBDQxcUVMpkMTU3Wl+z/UFVVJYDvX6heLTQ0DABQXV2FiRON329XVzez/flSqRQGg6H79oz7zX84UV4sFmPcuJA+92WLXntzI59riZub2zXbO3re1ff09LymTiyWmAVDPV/b1Xoe36qqKtTUGB/zv//9b/j73/9mVgsAtbU18PX1AwB4eHharLlafx/rvoqIMJ5MUFRUiMmTp8JgMODdd9+Gs7MLHn30f/p1W0RERERkW+1qHTZfqsaXZytQ3aY2ux7gqsDdNwVh5UR/OCvG/Mvi67KrR6C+vh6PPfYY7rzzTkRHRyMjI8N0Bn0PBwcHqFQqdHZ2Wr1may8sjMZTs8Oh1Y2caEAmFcFJPvCnj0gkQnx8ArKzM6HVaq95x/5q//znWhQXF+Hpp5/tdXWCwWC8dvXt9Da0z9iD8f9rNOa/KAwGAZIb2Eo02L325kY+19KKDYmVL7xvIZB5jSAYum9XbAorfvrTxxEfn2DxFkJCwkz/uy9fW38f677qOZmgZ8XA1q2bkJ+fi5///LkbXlFCREREREOjqlWFL89WYPOlanRozP/Wjfd3wb2TgzEv2htSMU8962E3wUBBQQEeeeQRzJo1C7/5zW8AWH6h39XVBUdHx16vDQcnuRQYYydjzJkzH+fPn8WePbuwdOkKs+tqtRpbt25CV1cXXF3d4O9v3BNeUlKMWbOurS0tLQEA0zvNfREYGNz9uaXXfFwQBFRWVlh897uvBrvXGyUWi01H9fXQ6XRoaWkxWzFxI3qm/l+trKwMgHHrQU8Q4ejoiClTpl1Tl5V1Ba2trVAoFP26z6F6rD08PODp6YWiogJ0dLTj44//H8LDI7B69e39vi0iIiIiGlpXqlrxnzMV2JdbB/0P3jcSAZgb7Y17bw7CpEBXHoNuwZifMQAAp0+fxj333IO7774br732mumJEBERgeLiYlOdXq9HaWkpIiIier1GgyMt7TYEBATi739/F4WF+ddcMxgMeOedP6GurhY/+tGPIZfLMWFCLLy8vPHttxvQ0dFuqu3oaMe3366Hl5c3YmJi+3z/48fHICAgEJs2bbgmBNq7d/d1tzj0vKtubbn+YPd6o7y8vFBaWgK1+vuv88iRQxZXS9yInJws5OZmm/67sbEBu3fvQGJiMtzc3DFhQhy8vLyxfv1X6OzsNNV1dLTj5Zd/jddff83qigVrhvKxjoyMQlFRAT799GM0NTXiF7/4Va+naBARERGR7egNAvbl1eOnX5zHA/89j4yca0MBB5kYdyUH4puHp+BPK+OQGOTGUMCKMf8XblVVFR5//HH87//+L+64445rrk2bNg3Nzc3YuHEjUlNTsXbtWgQGBiI6OhohISFWr9HgkMvleP31P+PZZ5/CT3/6E6SkLMGECbFoa2vD/v17kJubg9mz5+FHP7ofgHHf/DPP/Aovv/xrPPzwj7FiRRpEIhG2bduMhoZ6/P73b/ZrWb1IJMIvfvErvPjic3jssQexfPlK1NfXYuPGr+Hq6tbr5/YsJf/ii89wyy23YubMOddcH+xeb9TChYvxzjt/xi9/+TRSUpaivLwMW7Z8C3//gEG9H1dXVzz77JO46657IZFI8e2366HT6fDUU88CuPZxeeih+5Camga5XIGtW79FdXUVXn759/1+4T2Uj3VERBROnTqB9eu/wNy583HzzVMGdDtERERENHg6NDpsvVyDL89WoKLFfKu3j7McdyUHYfUkf7gq+7+l1B6NimDgpZdeQklJCT777DOza2VlZXjzzTdx8uRJAMDcuXPxwgsvmAalbdiwAW1tbXj99dfx+uuvmz7v5ZdfxurVq/GPf/wDr7zyCv7whz9gwoQJePfddyESiaBUKq1eo8ETHR2DTz/9D77++gscO3YY+/btgcGgR0REFF544SUsX552zWM+d+4CvP32+1i37mOsW/cxpFIp4uIm4oUXfovExOR+3/+MGbPwpz/9FZ988g/84x/vw8fHFy+88Ft88836Xj9vwYLFOHBgH3bs2Ipz586YBQND0euNWL36DrS2tmLr1k3461//jMjIaLz++p/x5ZefX/PO/Y2aNu1WTJgQhy+//BzNzc2Ii4vH7373BiZM+P4d+57H5d///ifWrfsEYrEIERGR+L//e9t0RGJ/DdVj3XMygUQixZNPPjPg2yEiIiKiG1fdqsLX5yrx7aUqtKvN5weM93HCvZODsSjGBzKJXSyOHzQiYaDnzdnI+vXr8dJLL2Hq1KlmwUBTUxPWrFkDjUaD+++/H3q9Hp988gmCgoKwfv16yOUjZ1N+XV3bcLdAREREREQ06mRWt+G/Z8qxJ8d8fgAAzIrwxL2Tg3FTMLcK/JCPT9+OXxyxKwb0ej0+/PBDvP/++1Zr1q1bh+rqamzduhWRkZEAgMTERDz44IPYtGkT7rzzTlu1S0RERERERINEbxBwqKABX5wpx7mKVrPrCqkYK+L9cM9NQQj1HJ4B8WPJiAwG1Go17rjjDuTk5GDVqlU4duyYxbrt27dj6tSpplAAAG699VaEh4dj+/btDAaIiIiIiIhGkevND/BykuOu5ECsnhQAdwfODxgsIzYYaG9vxzvvvINly5Zh/vz5ZjUtLS0oKyvD4sWLza7Fx8fjwIEDNuiUiIiIiIiIblR1qwpfnavEJivzA6J9nHDvzcb5AXIp5wcMthEZDDg7O2P37t29TievqakBAPj5mZ9P7uPjg/b2drS1tcHFpW97KoiIiIiIiMi2rlS14r9nKrA31/L8gJkRnrj35mDcPI7zA4bSiAwGxGLxdY8X6+joAAA4ODiYXVMoFACAzs5OBgNEREREREQjiM4g4GB+Pf57pgIXK63PD7j7piCEcX6ATYzIYKAvDAbDdWtseU48ERERERERWdeu1mHL5Wp8dbYCla1qs+veTnLcyfkBw2LUBgNOTk4AjPMIfqjnYz01RERERERENDwqWrrw1dlKbLlcjQ6N+fyA8T5OuHeycX6ATMI3d4fDqA0GAgMDAQB1dXVm12pra+Hq6gpHRy47ISIiIiIisjVBEHCx0jg/4EB+PQw/mB8gAjAr0gs/ujkINwVzfsBwG7XBgKurK4KDg3HlyhWza5mZmZg4ceIwdEVERERERGS/dHoD9ubW479nK5BZ3WZ23UEmRmq8P+66KQghHubz4mh4jNpgAABSUlLw73//GwUFBYiMjAQAHD16FEVFRXj44YeHuTsiIiIiIiL70NKlxaZL1fj6XAVq2zVm132d5bgrOQirJvnDVcn5ASPNqA4GHnnkEWzevBkPPPAAHnroIajVanz88ceIi4tDWlracLdHREREREQ0phU3duLLsxXYfqUGKp35gPg4fxfce3MQ5kd7Q8r5ASPWqA4GPD098fnnn+ONN97Au+++C6VSiQULFuBXv/oV5HL5cLdHREREREQ05giCgFOlzfjibAUOFzaaXReLgHnR3rjnpiBMCnTl/IBRQCQIgnD9MrpRdXXm+2uIiIiIiIhGC7XOgF1ZtfjibAXy6zvMrjvJJViVEIA7kwMR6KYchg7ph3x8XPpUN6pXDBCNVnv3ZmDDhi+Rn58LNzd3bNiwdbhbIiIiIiKyqKFDgw3nK7HxQhWaurRm14PclLj7piCkTvSDk5wvMUcjfteIhoGrqwtuv/0u1NXVYsOGr4a7HSIiIiIiM7m17fjv2Qrszq6FVm++0Dw52A0/uikIsyK9IBFzu8BoxmCAhtWTTz4KAHj//bV9+vhYMWXKLQCA/fv3DOrt2uvjSURERESDwyAI+K6gEV+eLcfpshaz61KxCIsn+ODum4Iwwa9vy9Rp5ONYSBpVOjraMWvWFDz99ONDej8ff/z/MHPmZFy4cM5qjcFgwIMP/gjz59+KioryMdkDEREREdmHDo0OX52twJp/nsJzm6+YhQLuDjI8fEsItj4yFa8uncBQYIzhigEaVXJysiEIAsaPnzCk9xMZGQUAKCwsQGJissWa7du3IC8vF/fd9wCCgoIBABqNBjqdzurtKhQKSCSSIe2BiIiIiKivKlq68PW5Smy+VI0Ojd7seoSXI350cxAWT/CFUta3v2Np9GEwQKNKbm42ACAmZqiDgWgAQFFRgcXrnZ0d+OijD+Ht7YP773/I9PE//vEV7N2bYfV23333/+GmmyYPaQ9ERERERL0RBAHnK1rxxdkKHMyvh8HCOXUzwj1xz81BmBrizuMG7QCDARpVcnJsEwwEB4+DQqFAYaHlF+Xr1n2CxsYGvPzy7+Ho6Gj6+GuvvYHXXntjWHsgIiIiIrJEqzcgI6cOX5ypQHZtu9l1pVSMFfF+uOumIIR58u9Le8JggEaV3NxsODo6Ydy40CG9H7FYjPDwSBQXF5pdq6gox4YNXyIhIREpKUsHdPt6vR46nQ46nQ6CIECtVkMkEkEul9usByIiIiKyD42dGnxzoQobLlShoUNjdt3PRYG7kgORluAPV6VsGDqk4cZgYJQQadoBvfkP8bCRyCHInW16l52dnSgrK0VCQqJNljNFREQiOzsTjY0N8PT0Mn38gw/+Bp1Oh1/84lcDvu1du3bg9ddfM/33ggUz4O8fgA0bttqsByIiIiIa23Jr2/Hl2Qrsyq6FxsJxgwkBrrjn5iDMi/aGlMcN2jUGA6OA88EXobz8GUSwsPlnmAgQQTXxx2if87rN7jMvLxcGg6HXbQRarRb//vc/+3yby5alIiAg0OK1qCjjHv/CwgLTi/KzZ0/j0KH9WLly9Q1tZ1i2LBXLlqVet24oeyAiIiKisUdvEHC4sAFfnK3AGQvHDUrEIiyI9sY9NwdhYoDrMHRIIxGDgZFOEOBw+d/D3YUZEYx9tc/+I2CjYSQ5OVkA0OuJBFqtFp9++lGfbzM5+WarwUBEhPFUgKKiQkyePBUGgwHvvvs2nJ1d8Oij/9OPzgduJPRARERERCNfu1qHLZer8fW5SlS0qMyuuymlWD0pALcnBcLPRTEMHdJIxmBgpBOJ0DXx/hG6YuB+m4UCwNUnEsRarXF0dMThw6cH5f56TgXoGf63desm5Ofn4uc/fw7u7u6Dch+joQciIiIiGrnKmrrw1bkKbLtSY/G4wXAvR9xzUxCWxvK4QbKOwcAo0D7ndXRMf9HuZwzk5GRBqVQiJGRoBw/28PDwgKenF4qKCtDR0Y6PP/5/CA+PwOrVt9vk/kdKD0REREQ0sgiCgFOlzfjybAUOFzZafPtwZoQn7k4OwtRQHjdI18dgYJSw9YvwkUatVqG0tARxcfGQSGyXdEZGRiEz8zI+/fRjNDU14tVX/wip1LY/NkPZQ0dHO15++UVUVpZj+fKVmD9/EX772xeg02nxk5/8FPPnLxyU+yEiIiKiG6fS6rEzqxZfnqtAQX2n2XUHmRgr4v1xV3IgQnncIPUDgwEaFfLycqHX63udLzAUIiKicOrUCaxf/wXmzp2Pm2+eYtP7H+oeNm/+Fvfeez8mTUrC//7vs9i7dzeef/43GD9+An772xcwZ848mwYxRERERGSutk2NDRcq8c2FKrSodGbXA10VuCM5CGkT/eGi5Es86j8+a2hUyMkxzhewdTDQcyqARCLFk08+Y9P7tkUPlZUVWLPmTkilUtx77/346KO/IzY2HgAwblwIWltb4OHhOaj3SURERER9c6myFV+ercDevHroDeYbBpKD3XDPTUGYHekFCY8bpBvAYIBGhTVr7sSaNXfa/H6XLl2BpUtX2Px+bdXD+PExyMhIx6xZc/Gf//wbMpkcW7duwpw581FQkAd3d48huV8iIiIiskyrN2Bvbj2+PFuBK9VtZtdlEhFSJvjinuQgxPjZ93ZjGjwMBojs2LJlqfjww/ewfftW3H773Zg+/Va89dab2LFjKx5//EkOqiEiIiKykcZODb69WIUN56tQ32E+dNzTUYbbEwNxW2IAvJzkw9AhjWUMBojsmFQqxVNPXbs94be//d0wdUNERERkf3Jq2/Hl2Qrszq6FRm++XSDWzxl33xSERTE+kEnEw9Ah2QMGA0RERERERDakMwg4lF+PL89V4lx5i9l1iQiYF+2Du28KxKRAV67ipCEnEgTB0rGXNMjq6sz3BxERERERkf1oVWmx+VI1vj5Xieo2tdl1N6UUqyYF4PbEAPi7KoehQxprfHxc+lTHFQNERERERERDqLChA1+drcSOzBqodAaz65HejrgrOQhLY32hlPGoaLI9BgNERERERESDzCAIOFLYiC/PVuBkabPZdRGA2ZFeuOumQEwe587tAjSsGAwQERERERENkna1DlsuV2P9+UqUN6vMrjvJJUhL8McdSYEIdncYhg6JzDEYICIiIiIiukHFjZ1Yf64S267UoFOrN7se4uGAu5KDsCLeD45ybhegkYXBABERERER0QAYBAHHipvw1dkKHCtuslhza7gH7kwOwvQwD4i5XYBGKAYDRERERERE/dCh0WH7lRp8da4SpU1dZtcdZRKsiPfDHcmBCPN0HIYOifqHwQAREREREVEflDZ14etzFdh2pQYdGvPtAsHuStyZHITUeD84K/hSi0YPPluJiIiIiIisMAgCTpQ04auzlThS1GixZlqoO+5KDsKt4Z6QiLldgEYfBgNEREREREQ/YNwuUIuvz1WgxMJ2AaVUjOXxfrgzORARXk7D0CHR4GEwQERERERE1K2sqQtfn6/E1svVFrcLBLopcWdSIFIn+sFVKRuGDokGH4MBIiIiIiKya33ZLjAlxLhdYGYEtwvQ2MNggIiIiIiI7NL1ThdQSMVYHmc8XSDKm9sFaOxiMEBERERERHbleqcLBLoqcHtSIFZO9IebA7cL0NjHYICIiIiIiMY8gyDgWHETvjpbgWPFTRZruF2A7BWDASIiIiIiGrPa1TpsvVKDDectbxfoOV3gjqRARHK7ANkpBgNERERERDTmFDV04utzFdieWYMurcHsOk8XIPoegwEiIiIiIhoT9AYBhwsb8fW5CpwsbbZYMy3UHXcmB2FGOLcLEPVgMEBERERERKNaS5cWWy5XY8OFKlS2qMyuO8okWB7vhzuTAhHm5TgMHRKNbAwGiIiIiIhoVMqra8dX5yqxM6sWap35doFx7krcmRyEFfF+cFbwpQ+RNfzpICIiIiKiUUOnN+BgQQO+OleJc+UtFmtuDffAnclBmB7mAbGI2wWIrofBABERERERjXiNnRpsuliNjRcqUduuMbvuJJdg5UR/3J4UiBAPh2HokGj0YjBAREREREQj1pXqNnx9rgIZOXXQ6gWz6+FejrgzKRDL4vzgKJcMQ4dEox+DASIiIiIiGlE0OgP25Nbh63OVuFLdZnZdLAJmR3rhzuRATB7nDhG3CxDdEAYDREREREQ0ItS0qfHNhUp8e7EaTV1as+tuSinSEozbBQJclcPQIdHYxGCAiIiIiIiGjSAIOFfRgq/PVeJAXj0s7BZAjK8z7kwOREqMD5QybhcgGmwMBoiIiIiIyOa6tHqkZ9bg6/OVKKjvNLsuEYuwcLw37kgKxKRAV24XIBpCDAaIiIiIiMhmSpu6sOF8JbZeqUa7Wm923ctJjjWTArB6kj+8nRXD0CGR/WEwQEREREREQ8ogCDha1Iivz1XiWHGTxZrEQFfcmRyIedHekEnENu6QyL4xGCAiIiIioiHR0qXF1is12HC+EhUtKrPrCqkYSyb44o6kQMT4OQ9Dh0QEMBggIiIiIqJBllPbjvXnKrEzuxZqncHseqCbErcnBmDlRH+4OciGoUMiuhqDASIiIiIiumFavQH78+rx9blKXKhstVhzS5gH7kwKxK3hnpCIOUyQaKRgMEBERERERANW167GNxeq8O2lajR0aMyuOyskWBHvj9sTAxDq6TgMHRLR9TAYICIiIiKifhEEAecqWrD+XCX25zdAbxDMaqK8nXBHUgCWxPrBUS4Zhi6JqK8YDBARERERUZ90avRIz6rB+vOVKKjvNLsuEQHzor1xR3IgkoPcIBJxuwDRaMBggIiIiIiIelXc0IkNFyqx7UoNOjR6s+uejjLcNikAqycFwNdFMQwdEtGNYDBARERERERmdAYBhwsasP58JU6WNlusSQpyxR1JgZgX7Q2ZRGzbBolo0DAYICIiIiIik8ZODTZfqsbGC1WoaVObXVdKxVgS64s7kgIx3td5GDokosHGYICIiIiIyM4JgoBLVW1Yf74Se3ProNWbDxMM8XDAmsQApMb7w0XJlxFkf8QdNVDkb4Ws/Ah03nHonPIsIB4bgzX5E01EREREZKdUWj12Zddi/fkq5NS2m10XAZgZ4Yk7kgMxLdQDYg4TJDsjUrdAUZAORd4myCqOQiQYAACK4gxox82GNnDaMHc4OBgMEBERERHZmbKmLtMwwVaVzuy6m1KKtIQArEkMQKCbchg6JBpGui7Ii/dCmbcJ8uJ9EBk0ZiWC1AF656BhaG5oMBggIiIiIrIDeoOAI0WNWH++EseLmyzWxPu74I6kQCyM8YFCymGCZEcMOsjKjxjDgIJ0iLXmK2gAQO/oC3X0SqjifwyDa7CNmxw6DAaIiIiIiMawpu5hgt9crEJVq/kwQYVUjEUxPrgjKRBx/i7D0CHRMBEESGvOQpn7LRT52yDuqrdYZpC7Qh25FOro1dAGTR8zcwWuxmCAiIiIiGiM6RkmuOF8JfZYGSYY6KbE7YkBSJ3oD3cH2TB0STQ8JA3ZUOZugiJvMyRtZRZrBIkCmrCFUI1fBU3IPEA6trfUMBggIiIiIhojVFo9dmbVYsMF68MEZ0R44vbEQEwP5zBBsh/i1jIo8jZDmbcJ0oZsizWCSAxt8CyootOgiVgCQeFq4y6HD4MBIiIiIqJRrrixE99cqMLWK9VoV+vNrrsppVg50R+3JQYg2N1hGDoksj1RZx0U+dugzNsEWfUZq3Va/5uhil4FddQKCI4+Nuxw5GAwQEREREQ0CukMAr4raMCG85U4WdpssSbe3wW3JwVg4XgfKGVjb1800Q+J1K2QF+6EMm8zZOXfmY4X/CGdZwzU0augil4Jg1uojbsceRgMEBERERGNIvXtamy6VI1vL1ahtt38GDWFVIyUGB+sSQpEPIcJkj3QdUFesg/K3E2Ql+yDSG8+ZBMA9C7BUEenQTV+FfResTZucmRjMEBERERENMIJgoCz5S3YcL4K+/ProTeYDxMc567EmsRArIj3gxuHCdJYZ9BBVvYdlHmbIS/cafV4QYODF9RRqVBFr4LO/2aAczUsYjBARERERDRCtat12JFZgw0XqlDU0Gl2XSwCZkV44fakAEwN5TBBGuMEA6RVp6HM2wxFwTaIuxoslhlkztBELoUqehW0wTMAMV/2Xg8fISIiIiKiESa3th0bL1QhPasGXVrzPdIeDjKsmuSP1ZMCEOA6to9RIzsnCJDWX4EibxMUeVsgaa+0XCZRQBM633i8YOh8QMohm/3BYICIiIiIaATQ6AzYl1ePDecrcaGy1WJNUpArbk8MxLxob8ilYht3SGQ7kuZCKHI3QZG/BdKmfIs1gkgC7biZUEWvgiZ8sV0dLzjYGAwQEREREQ2jyhYVvrlYhS2XqtHUpTW77iiTYGmcL25PDESUj9MwdEhkG+L2SijytkKRtxmyuotW67QBU6CKToM6cgUER28bdjh2MRggIiIiIrIxvUHAseJGbDhfhaNFjTAfJQhEejvi9sRALI3zhZOcf7bT2CTqaoSiYDsUeZsgqzwJkcWfBkDnFQfV+FVQR62EwTXYxl2OffwNQ0RERERkI42dGmy5VI1vLlahqtX8SDWpWIT50d64PSkQSUGuEHGYII1BIk075EU7ocjdDHn5dxAZdBbrdG5hUEenQR29CnrPaBt3aV8YDBARERERDSFBEHC+ohUbL1Rib249dBaOGvR3UeC2xACsnOgPLyf5MHRJNMR0KshL9hqPFyzeC5HePBgDAL2TP9RRK6EenwadzyQeL2gjDAaIiIiIiIaA8ajBWnxzsRIF9eZHDYoATA/3wO2Jgbg13BMSMV8A0Rij10JWftgYBhTuhFjbbrHMoPSAOnI51NFp0AZOA0QcrGlrDAaIiIiIiAZRTm07Nl6oxM6sWotHDboppUhLMB41GOzOI9VojBEMkFWdhCJvCxT52yBWNVosM8icoAlfDHV0GjTjZgMSmY0bpasxGCAiIiIiukFqnQF7cuqw8UIlLlW1WayZFOiKNYkBWDDeBwoeNUhjiSBAWnexOwzYAkl7leUyiQKa0PlQRadBE7oAkDEYGykYDBARERERDVBpUxe+uVCFbVeq0aIyH6DWc9TgbZMCMN7XeRg6JBo6ksZcKPI2Q5G3GdKWYos1gkgC7biZUEWvgiZ8MQSFq22bpD5hMEBERERE1A86g4BDBQ345kIlTpQ0W6yJ8nbCmsQALIn1hbOCf3LT2CFuKYEifyuUeZshbciyWqcJmAb1+DSoI5ZBcPS2YYc0EPwtRURERETUB7Vtamy6VIVNl6pR164xuy6TGI8avCMpEJMCedQgjR3ijmoo8rdBkbcZsppzVuu0PpOMxwtGpcLgEmjDDulGMRggIiIiIrLCIAg4WdKEjReq8F1BA/TmJw0i0E2JNZMCkDrRDx6OPGqQxgZRVyMUBTugyN8MWcVxiGDhyQ9A5xFtDAOiV0LvHmHjLmmwMBggIiIiIvqB5k4ttl6pxjcXq1DerDK7LhYBMyO8sCYxALeEeUDM1QE0BojUrZAX7YYybxNkZd9BJOgt1uldxkEdvRKq6DTovWIBPv9HPQYDREREREQABEHAhYpWbLxYhb25ddBaWB7g7SRHWoI/ViX4w99VOQxdEg0ybRcUxXugyN8Mecl+iPRqi2V6Rz+oo1ZAHZ0GnV8yw4AxhsEAEREREdm1drUOOzJr8c3FShTUd1qsmRLijtsTAzA70gtSCY8apFFOr4a89BAUeZugKMqASGf5eW9QuEMduRzq8WnQBkwDxBIbN0q2wmCAiIiIiOxSVk0bNl6owq6sWqh0BrPrbkopVsT7Y/Ukf4R6Og5Dh0SDyKCDrOKo8XjBwp0Qq1ssl8mcoYlYDHXUSmjGzQYkMhs3SsOBwQARERER2Y0urR67s2ux8UIVsmraLdZMCnTFmsQAzI/2hlLGd0hpFBMMkFWdgiJvCxQF2yDuarBcJlFAHbYI6uiV0ITOA6QONm6UhhuDASIiIiIa8/LrO/DthSpsz6xBh8Z8oJqjTIKlcb5YkxiAaB/nYeiQaJAIAqS1F4xhQP4WSDqqLZeJZdCEzDWGAWGLIMj5vLdnDAaIiIiIaExS6wzYm1uHby5U4UJlq8Wa8T5OWJMYgMWxvnCS809jGqUEAZKGLCjyt0KZtwWS1hLLZSIxtEEzoI5OhTpiKQSlh40bpZGKv/2IiIiIaEwpaezENxersP1KDVpUOrPrCqkYKTE+WJMYgDh/F4g4XZ1GKUlTART5W6DI2wJpU57VOk3ANKijV0IduQyCo48NO6TRgsEAEREREY16Wr0BB/Ib8M3FKpwubbZYE+7piNsSA7AszheuSg5Uo9FJ3FoKRf5WKPK2QFZ/xWqd1jcR6ug0qCNXwOASaMMOaTRiMEBEREREo1ZFSxc2XazGlsvVaOzUml2XSUSYH+2N2xIDkBzkxtUBNCqJ26ugKNgORd5myGrOWa3TeU2AOmolVNErYXALs12DNOoxGCAiIiKiUUVnEHCksAEbL1TheHETBAs1we5K3DYpACvi/eDhKLd5j0Q3StRZ3x0GbIGs6iREFp/pgM49AuqolVBHr4Tec7yNu6SxgsEAEREREY0K1a0qbL5kXB1Q264xuy4RAXOijKsDpoS4Q8zVATTKiFRNUBSmQ5G3FbKKIxAJBot1epdg48yAqJXQeccDfK7TDWIwQEREREQjlt4g4FhxIzZeqMLRokYYLLxp6u+iwOpJAVg50Q/ezgrbN0l0A0SaNsiLdkGRtxXysoMQGcwHZgKA3tEP6qjlUEenQed3E8MAGlQMBoiIiIhoxKlrV2PzpWpsvlSN6ja12XWxCJgR7ok1iYG4JcwDEjFfJNEoou2EongPFPlbIC/ZD5He/DkOAAalJ9RRK6COSoU2YCoglti4UbIXDAaIiIiIaEQwCAKOFzfh24tV+K6gAXoLqwN8nOVIm+iPtAR/+Lsqbd8k0UDpVJCX7ocibwsUxXsg0nVZLDMo3KCOWAJ1dBq0QbcCYr5ko6HHZxkRERERDav6djW2XqnBpotVqGw1f+dUBGB6uAdumxSAGRFekHJ1AI0Weg3kZYegyNsCedFuiLXtFssMMidowlOgjk6DZtxsQMKBmWRbDAaIiIiIyOYMgoCTJU345mI1DhU0QG9heICXkxxpE/2QlhCAQDeuDqBRQq+FrOIIFHlboSjaCbG6xWKZIHWAOmwh1FGp0ITOA6QONm6U6HsMBoiIiIjIZho6NNhyuRqbLlWjskVlseaWUA+snuSP2ZFekErENu6QaAAMesgqjxu3CRTugFjVZLFMkCigCZ1nPF4wbCEgc7Rxo0SWMRggIiIioiFlEAScKmnGt5eqcCDf8uoAT0cZUif6Y1WCP4Ld+c4pjQKCAbKqU1Dkb4EifwfEXXWWy8QyaMbNhjo6FZrwxRDkLjZulOj6GAwQERER0ZBo6NBga/fqgAorqwOmhLjjtkkBmBPlBRlXB9BIJwiQ1pyFIn8rFPnbIOmotlwmkkAbPBPqqFSoI5ZAULrbtk+ifmIwQERERESDpi+rAzwcZEid6IdVCQEY58HVATTCCQKkdReN2wTyt0HSXmG5DCJog6YbtwlELoXg4GXjRokGjsEAEREREd2w+g4Ntl1ndcDkntUBkV6QS7k6gEYwQYCkIQvKvC1Q5G+FpLXEaqk2YCpUUalQRy6H4ORrwyaJBg+DASIiIiIakJ6TBb69WI2DVk4WcHeQITXeD6smBSCEqwNohJM0ZJu2CUibC6zWaf2SjSsDopbD4Bxoww6JhgaDASIiIiLql/qrZgdYO1lgcog7Vif4Y26UN1cH0Igmaco3bROQNuVardP6TDLODIhaAYPrOBt2SDT0GAwQERER0XUZBAEnulcHHOpldcCKeD+sSvBHqCePYaORS9xcBGX+Vijyt0LakGW1TucVC3XUSqiiVsDgHm7DDolsi8EAEREREVlV167G1ss12HypCpWtaos1U0LcsZqzA2iEE7eUmLYJyOovW63TecZAHbUC6qhU6D2ibNgh0fBhMEBERERE19AbBBwvbsK3F6twuLABevPFAfB0lGFFvD9WJfjzZAEascSt5VAUbIMifytktRes1uncI7u3CaRC7xVjww6JRgYGA0REREQEAKhpU2PL5WpsvlSNmjbLqwOmhRpXB8yO9IJMwtUBNPKI2yqhKNgORf4WyGrOWa3Tu4ZCFb2yOwyIBUQiG3ZJNLIwGCAiIiKyYzqDgKNFjfj2YhWOFjXCwugAeDrKsHKiP9IS/BHsztUBNPKI26u6w4BtkFWftlqndxkHdbRxZYDOeyLDAKJuDAaIiIiI7FBliwqbL1dj6+Vq1LVrzK6LANwS5oFVkwIwO8ITUq4OoBFG3FEDecF2KPO3QVZ10mqd3jnINDNA55vIMIDIAgYDRERERHZCpzfgUEEDvr1UjRPFTbCwOAA+znKkTvRH2kR/BLopbd4jUW9EnXVQFOwwzgyoPAGRxWcxoHfyNx0tqPO7iWEA0XUwGCAiIiIa48qaurDpUjW2XalGY6fW7LpYBNwa7olVCQGYEeEJqZgvomjkEHXWQ1GY3h0GHIdIMFis0zv6QR21HOqoldD53wSIuMqFqK8YDBARERGNQRqdAQfy6/HtpWqcLm22WOPnokDaRH+kTvSDvytXB9DIIepq6A4DtkFWcbSXMMAXmshlUEelQhswhWEA0QAxGCAiIiIaQwobOrDpYjV2ZNagRaUzuy4RAbMivbBqUgBuCfWAhKsDaIQQqZq+DwPKj0Ak6C3WGRy8oI5c3h0GTAXEEht3SjT2MBggIiIiGuVUWj0ycuqw6VI1Lla2WqwJdFNiVYI/UuP94O2ssHGHRJYZw4Cd3WHA4d7DgIilxjAgcBog5ssYosHEnygiIiKiUSqnph3fXqrCzqxadGjMX1BJxSLMjfLGqgR/TAl1h5gD2GgEMIYBu6Ao2GpcGWAwX9kCAAalB9QR3dsEgm5hGEA0hPjTRURERDSKtKt12J1di02XqpFV026xJtTDAasmBWB5nC88HOU27pDInEjVDHnRLijztxpXBlgLAxTuUEf2rAyYDkhkNu6UyD4xGCAiIiIa4QRBwKWqNmy6WIWMnDqodOaD2BRSMRaM98aqhAAkBblCxNUBNMx6wgBF/jbIy7/rJQxw694msALaoBkMA4iGAYMBIiIiohGquUuLHZk12HypGoUNnRZron2csCrBH0tifeGq5AsqGl7GMGA3FPlbIS8/DJHB/HhMoCcMWAJ15Apog2cyDCAaZgwGiIiIiEYQgyDgTFkzNl2sxv78emj1glmNo0yClAk+WDUpAHF+zlwdQMPKFAYUbIO87Lvew4DwJcaVAcEzAAm3uRCNFAwGiIiIiEaAunY1tl0xrg6oaFFZrJkY4IJVCf5YFOMLRzmPaKPh058wQBO+GOqoFdAEz2QYQDRCMRggIiIiGiY6g4CjRY3YdLEKR4saYWFxAFyVUiyN9cWqhABE+TjZvkmibtfODOh9m4AmfDHUkcuhGTeLYQDRKMBggIiIiMjGypu7sOVyNbZerkF9h8ZizeRxbliZ4I/50T5QSMU27pDIqM8DBOWu0EQsYRhANErdUDDQ2dmJPXv2ICsrC+3t7fj973+PS5cuQSKRIC4ubrB6JCIiIhr11DoDDuTVY9PlapwubbZY4+UkR2q8H1ZO9Mc4DwfbNkjUTaRqgqJwFxQF23o/WlDuCk3EYqgjVzAMIBrlRIIgWFi0dn179uzBb37zG7S2tkIQBIhEImRlZeGdd97B2rVr8eCDD+L5558f7H5Hrbq6tuFugYiIiIZBXl07Nl+qRnpWLVpV5i+wxCLg1nBPrErwx4xwT0glXB1AtmcMA3Z2hwFHGAYQjRE+Pi59qhvQioGzZ8/iF7/4BXx9ffGzn/0MmZmZ2LJlCwBgxowZ2L17Nz799FPExsYiNTV1IHdBRERENGq1q3XYnVOHzZeqkVlt+c2BQFcF0hICsCLeD74uCht3SASIuhqhKNoJRf5248oAQW+xjmEA0dg3oGDggw8+gJubGzZs2ABPT0+8//77pmtTp07Fl19+idTUVPznP/9hMEBERER2QRAEXKxsxeZL1cjIqYNKZzCrkUlEmBfljbQEf0wOcYeYxwySjYm6GqAoTDeGARVHrYcBHCBIZFcGFAxcuHABaWlp8PT0tHjdzc0NKSkpplUERERERGNVY6cGOzJrsflSFYobuyzWRHo7Ii0hAEtjfeHuILNxh2TvRJ31328TqDjWaxigjlgCTeRyHi1IZGcGFAzo9XqIrpNw63Q6aLWWjzAhIiIiGs30BgHHS5qw5VI1DhY0QG8wH9nkIBMjZYIvViX4I97f5bp/OxENJlFnnTEMyN8GWeUxiATzFSwAYFC4Qx2xBOqoFdAGzQAkDK6I7NGAgoHx48fju+++g1arhUxm/sujo6MDhw4dQkxMzA03SERERDRSVLaouo8ZrEZtu+VjBhMCXJGW4IeFMT5wkvNkaLIdUUetcZtAwTbIKk9YDwOUHsYwIHIFtEG3MgwgooEFAz/+8Y/x3HPP4Wc/+xmee+65a1YGZGVl4fXXX0dVVRV+/vOfD1qjRERERMNBrTPgYH49Nl2qxikrxwy6O8iwLM4XaQn+iPBysm2DZNfEHdWQF+yAomA7ZJUnIYLlA8cMSk+oI5YaVwYE3sIwgIiuMeDjCt9++22sXbv2mmVxUqkUOp0OgiDgxz/+MX7zm98MWqODIT09HZ999hn++9//mj528eJFvPLKKygqKkJMTAzeeOMNREREXPdaf/G4QiIiotElt7YdWy5bP2ZQBGBamAdWJfhjdqQXZDxmkGxE3F4JRU8YUHXKap3BwcsYBkSugDboFkDMFSxE9qavxxUOOBgAgPPnz+Obb75BZmYm2tra4OjoiJiYGKxevRrTpk0b6M0OOr1ej3/96194++23kZCQgC+++AIAoFarsWDBAjz33HNYtmwZ1q5diwMHDmDDhg29XhsIBgNEREQjX5tKh13ZtdhyuRpZNe0Wa/xdFFg50R+pE/3g76q0cYdkr8St5VAUdocB1Wes1hkcfKCOXAZ15DLjygCxxIZdEtFI09dg4IZiw6SkJCQlJd3ITdjEX/7yF1y8eBGPPvoojh07Zvr48ePH4ezsjFWrVgEAnnjiCfzrX/9CTk4OqqurrV7j7AQiIqKxQxAEnC1vweZL1diXVw+1lWMG50Z5I22iP6aE8phBsg1xS0n3yoBtkNVesFqnd/SDJrJ7m4D/FIYBRNRvdrGe6MEHH4Svry+++eabaz5eVFR0zdYAiUSCcePGoaCgALW1tVavMRggIiIa/Wrb1NieWYMtl6tR3qyyWBPl7YSVCf48ZpBsRtxcBEXBduPKgLpLVuv0Tv5QRy6HOnI5dAGTARG3shDRwA0oGIiNje1TnUgkQmZm5kDuYlD5+vpa/HhnZyccHByu+ZiDgwNUKlWv14iIiGh00uoN+K6wEVsuVeNYcSMsnDIIJ7kEiyf4YmWCP+L8nHnMIA05SVO+MQzI3w5pg/W/nfXOQcYwIGo5dH7JDAOIaNAMOBiw9I+kSqVCZWUlurq6EB8f3+cAYbhYeqHf1dUFR0fHXq8RERHR6FLY0IHNl6qRnlmLpi6txZrkYDekTfTHgvHeUMq4FJuGkCBA0phrWhkgbcyxWqp3GQd1VPfKAN8kgEEVEQ2BAQUDP1ySfzW9Xo9//vOf+Pvf/45XX311oH3ZRERExDXDBPV6PUpLSxEREQEHBwer14iIiGjka1frkJFTh62Xq3GpyvIQYG8nOVbE+yF1oj9CPBws1hANCkGApCHr+zCgKd9qqd41tDsMWAGdTwLDACIacoM+Y0AikeCRRx7B8ePH8fbbb+PTTz8d7LsYNNOmTUNzczM2btyI1NRUrF27FoGBgYiOjkZISIjVa0RERDQyCYKAcxUt2HK5Bntz6qCyMEhQIhZhVoQnVk70x/RwT0jFfNFFQ0QQIK27BEXBdsgLtkPaUmy1VOceAXXkCqgjl0PvHccwgIhsasiGD8bFxeHzzz8flNt66aWXUFJSgs8++8zsWllZGd58802cPHkSADB37ly88MIL8PT0vO7tKpVK/OMf/8Arr7yCP/zhD5gwYQLeffddiESiXq8RERHRyNIzSHDr5WqUWRkkGObpgJUT/bEszg9eTnIbd0h2QzBAWnOu+zSBHZC0lVkt1XnGdB8tuBx6zxiGAUQ0bESCIFgYu3NjBEHAvffei8LCQhw/fvyGbmv9+vV46aWXMHXqVLNgoKmpCWvWrIFGo8H9998PvV6PTz75BEFBQVi/fj3k8pHzj35dneUljERERDQwWr0B3xU0YMvlGquDBB1lEiya4IOVE/2REODCgJ+GhmCAtOq0cZtA4Q5I2qusluq8YqGOWgF1xDLoPbkSlYiGlo+PS5/qBrRi4D//+Y/FjwuCgI6ODhw6dAjnzp1DWlraQG4egHFP/4cffoj333/fas26detQXV2NrVu3IjIyEgCQmJiIBx98EJs2bcKdd9454PsnIiKikSm/rgNbLlcjPasWzdYGCQa5InWiPxbG+MCBgwRpKBh0kFWegKJgB+SFOyHprLFaqvWZBHXkMmgil0HvznlVRDTyDCgY+P3vfw+RSITeFhvEx8fjueeeG1BTarUad9xxB3JycrBq1SocO3bMYt327dsxdepUUygAALfeeivCw8Oxfft2BgNERERjRJtKh13ZtdhyuRpZNe0Wa3oGCa6I90OoJ08RoiGg10JWcdS4MqBoF8RdDVZLtX7JxqMFI5fB4BpiwyaJiPpvQMHAG2+8YfHjIpEIMpkMERERN3RUoVqtRnt7O9555x0sW7YM8+fPN6tpaWlBWVkZFi9ebHYtPj4eBw4cGPD9ExER0fAzCAJOlzZjy+VqHMhvgNrKIMHZkV5IjffjIEEaGno15GWHjQMEi3ZBrG6xWqr1n9y9TWApDC5BNmySiOjGDCgYWL169WD3cQ1nZ2fs3r0bUqn19mpqjMu1/Pz8zK75+Pigvb0dbW1tcHHp254KIiIiGhkqW1TYdqUa267UoKpVbbEmwssRaQn+WBrrCw/HkTNTiMYIbRfkZQegyN8OecleiDWWZ0UJIjG0gdOgjlwOTcQSGJz8bdsnEdEgGbJTCW6EWCyGWCzutaajowMA4OBgfuawQqEAAHR2djIYICIiGgVUWj3259djy+UanC5ttljjJJdgSawvUif6I87PmYMEaVCJNO2Ql+yDvGAHFCV7IdJ1WawTRBJog2cYTxMIXwLB0dvGnRIRDb4+BQPTp08f0I2LRCIcPXp0QJ97PQaD+XLCH7peuEBERETDRxAEZFa3YcvlGuzOqUW7Wm+xbvI4N6xM8Me8KG8oOUiQBpFI3QJ5cQYUBemQlx6ASG95hYoglkEzbpZxZUB4CgSlh407JSIaWn0KBiy9Kz/cnJycABjnEfxQz8d6aoiIiGjkaOjQID2rFlsvV6OwodNijb+LwjhIcKIfgtxG3t8hNHqJuhqhKNoFecEOyMsPQ2SwfLKFIFFAEzLXGAaELYSgcLVxp0REttOnYGDfvn1D3Ue/BQYGAgDq6urMrtXW1sLV1RWOjpxITERENBLo9AYcKWrElss1OFLUCL3B/GQjuUSEedHeSJ3ojykh7hBzqwANElFHLRRFO6HI3w5Z5XGIBMurUwSpI9RhC6CJWAZ16HxAzjeZiMg+jMgZA33h6uqK4OBgXLlyxexaZmYmJk6cOAxdERER0dXy6zuw9XI1dmbVorHT8juzcf4uSI33Q8oEH7gqZTbukMYqcVslFIU7IC9Ih6zqJESwfMy2Qe4CTdhC48qAkDmAlCtUiMj+DDgYaG9vx7Fjx9DY2Ai9Xg9B+P6XrVarRXNzM44cOYL169cPSqOWpKSk4N///jcKCgoQGRkJADh69CiKiorw8MMPD9n9EhERkXWtKi12Z9dh65UaZFZbnubu4SDD0jhfpMb7I8qH78rS4BA3F0FRuAOKgh2Q1V6wWmdQuEMdsRiaiGXQjJsJSBQ27JKIaOQZUDBQUFCAH//4x2hqarrm44IgXDMhWC4f2uODHnnkEWzevBkPPPAAHnroIajVanz88ceIi4tDWlrakN43ERERfU9vEHCqtAlbL9fgQH49NHrzd2clImBGhBdS4/0wM8ITUgmHBNMNEgRIGnNNYYC0IctqqcHBB+qIJVBHLoM28BZAwtUpREQ9BhQMvPvuu2hsbMTy5csxbdo0vPfee4iPj8e8efOQn5+PDRs2wMnJCdu2bRvsfq/h6emJzz//HG+88QbeffddKJVKLFiwAL/61a+GPJQgIiIioLSpC9uvVGPblRrUtmss1oR7OWLlRH8sjfWFlxP/faYbJAiQ1l82HitYsAPS5gKrpXrnQKgjlkITuQxa/8mAmKdaEBFZIhKu3gPQRzNmzEBISAi++OILAMDzzz+P8vJy/Pe//wUAnD9/Hvfddx8ee+wxPPXUU4Pb8ShVV2d5KSUREdFo06HRYW9OPbZeqcb5ilaLNc4KCVJifLFyoh/i/F2uWVFI1G+CAdLqs1AU7ICiMB2StjKrpXrXUKgjl0EduQw63ySAzz0ismM+Pi59qhvQioGWlhYkJyeb/jsmJgYZGRmmrQRJSUmYM2cO9u7dy2CAiIhoDDAIAs6Vt2DrlRrsy61Dl9ZgViMCMDXUHanx/pgT5QWljO/O0g0w6CCrPAFFwQ7IC3dC0lljtVTnMR7qyKVQRy6H3iuWYQARUT8NKBhwdnaGTqcz/XdoaChUKhWKi4sRHh4OAIiMjMSpU6cGp0siIiIaFpUtKmy/UoNtmTWobFFZrAl2V2JFvB+Wx/nB31Vp4w5pTNGrIS87DHlhOhRFuyFWNVot1fokQBOxFOrIZdB7RNmwSSKisWdAwUB8fDwOHTqE5557DnK5HJGRkRAEAWfOnDEFA2VlZVw2SERENAp1afXYn1ePrZercbqsxWKNg0yMBeN9kDrRD8lBbvw3nwZO2wV56X7jyoCSvRBrrG+/1PrfDHXEMqgjl8LgGmLDJomIxrYBBQP33HMPnnzySaxZswa/+93vkJycjJiYGLz11lsQBAF1dXXIyMjAtGnTBrtfIiIiGgKCIOBCRSu2XanBntw6dGj0FuuSg92QGu+HBeN94CjnVgEaGJG6FfLiPVAUpkNeuh8ineXVKIJIDG3gNKgjlkETsQQG5wAbd0pEZB/6NHxwzpw5SEtLQ1paGiIjIwEA69atw3vvvYc//vGPWLJkCU6cOIHHHnsMarUagiDA1dUVn376KeLj44f8ixgNOHyQiIhGoupWFbZn1mD7lRqUNVt+cebvosDyeD+siPdDsLuDjTuksULU1QBF0W7IC3ZAXn4YIoPWYp0glkEbPMM4QDB8MQQHLxt3SkQ0dvR1+GCfgoGkpCSoVCqIRCIkJCRg1apVWL58ORwcHKDX6+HgYPwjobKyEhkZGVAqlZg7dy78/Pxu7KsYQxgMEBHRSKHS6rE/vx7bLtfgVGkzLP0hoJCKMT/aGyvi/TA5xB1ibhWgARC3V0JeuAuKwnTIKo9DJJgPrQQAQaqEJmSucWVA2AIICjcbd0pENDYNajDQ2dmJjIwMbNu2DceOHYNOp4NMJsO8efOwevVqzJ49GxIJlxP2hsEAERENJ0EQcLGyFVuv1GBPjvWtAomBrlge74dFMT5wVgxoxyHZOXFLMRQF6VAU7oCs5pzVOoPMGZqwhVBHLoUmZB4gc7Rhl0RE9mFQg4GrNTU1YceOHdi6dSsuXLgAAPDw8EBqaipWrVqF2NjY/ndrBxgMEBHRcKhuVWFHZi22Xam2ulXA11mOFfF+WBbnh1BPvjijfhIESBpzoChMh6JgB6QNWVZLDUoPqMNToIlYBs24mYBEYcNGiYjsz5AFA1erqKjA1q1bsX37duTl5UEkEiE6Ohq33XYbVqxYAW9v74He9JjDYICIiGxFpdVjX149tl2pweletgrMjfJCarw/Joe4QyLmVgHqB8EAae0F4/DAgnRIW4qsluqd/IzHCkYshTZwGiDmShQiIluxSTBwtezsbGzbtg27d+9GaWkppFIpZsyYgX/84x+DcfOjHoMBIiIaSn09VSAhwBWpE7lVgAbAoIOs6iTkBelQFO2EpL3KaqneNRTqSGMYoPNLBkRiGzZKREQ9bB4M9FCr1fjXv/6FDz/8ECqVCllZ1peT2RMGA0RENBQqW1TYkVmD7Zk1KO9lq8DyeD8s51YB6i+9GvKyw5AXpkNRtBtiVaPVUp1nDNQRS6GOXAa9VyzAgZVERMOur8HAoLxVoFKpsH//fqSnp+O7776DSqWCq6sr0tLSBuPmiYiI6CqdGj325dVh25UanClrsVijkIoxr+dUgXHcKkD9oOmAvHS/cZtA8V6Ite1WS7W+ScbhgRFLoXePsGGTREQ0mAYcDGg0Ghw8eBDp6enYv38/VCoVJBIJZs+ejVWrVmHevHmQyWSD2SsREZHdMggCzpQ1Y/uVGuzLq0eX1vKxb0lBrlgR74cF47lVgPpOpGqCvHgPFAXpkJcdhEivtlgniMTQBkw1zQwwuATauFMiIhoK/fqLQavV4rvvvkN6ejr27duHzs5OCIKAiRMnYtWqVVi+fDk8PDyGqlciIiK7U9rUhe2ZNdhxpQbVbZZfrPm7KLCse6tAiIeDjTuk0UrcXgV50S4oCndCVnEMIsHyXApBLIMmeCY0kUuhDkuB4Mjh0kREY02fgoFDhw4hPT0de/bsQXt7OwRBQGBgIO677z6kpaUhIoJLx4iIiAZLm0qHjNw6bL9Sg4uVrRZrlFIxFoz3xop4f9w0zg1i7uemPhA3FxmPFSzcCVnNWat1gtQBmtD5UEcshSZ0PgSFqw27JCIiW+vT8MEJEyYAABwdHbF48WKkpaXhlltuGfLmxhIOHyQiot7oDAJOlDRh+5UaHMyvh0Zv+Z/nm8e5YXmcH+aP94aTnFsF6DoEAZKGLCgKdkBRtBPShmyrpQaFGzThKVCHL4EmZDYg5eoTIqLRblCHD86YMQNpaWlISUmBUqm8ocaIiIjoewX1Hdh+pQbpWbWo79BYrAlyU2J5nB+WxfsiyI0v1ug6BAOk1WegKNwJRWE6JK2lVkv1jn7QRCyGOmIptIG3ABLOhyIiskeDflwhWcYVA0RE1KOpU4Nd2catAtm1lie+O8klWBjjgxVxfkgMcoWIWwWoN3oNZBVHu8OAXRB31VkvdQ2FOtI4PFDnlwyIxDZslIiIbMmmxxUSERFR7zQ6Aw4XNWL7lRocKWqE3mCey4tFwNRQD6yI88OcKC8oZZJh6JRGDW1n97GCO43HCmosz6MAAJ1XLNQRS6GOWAK9VyzAoImIiK7CYICIiGiICIKAzOo2bLtSg4ycOrSodBbrwr0csSLOD0vjfOHjrLBxlzSamI4VLNwJeekB68cKQgSd/81QRyyBOmIJDG5htm2UiIhGFQYDREREg6y6VYX0rFrsyKxBcWOXxRo3pRRLYn2xPN4PE3yduVWArBK3V0Je2H2sYOXxXo4VlEIbNAPqiCXQhKfA4ORn406JiGi0YjBAREQ0CDo1euzLq8P2zFqcKW2GpQE+UrEIMyM8sSLeD7eGe0Im4d5uskzSVAB5z7GCteet1glSJTQhc7uPFVwAQelusx6JiGjsYDBAREQ0QHqDgNNlzdiRWYN9ufVQ6QwW6+L8XbA8zhcpMb5wd+TUd7JAECCtuwh54U4oCndC2pRntdSgcIMmbJFxZcC4OYCMJ1UQEdGNYTBARETUT4UNHdiRWYv0zBrUtls+YtDXWY6lcX5YFueLCC8nG3dIo4JBB1nlCSgK0yEv2gVJe5XVUr2THzQRS6AOXwpt4DQeK0hERIOKwQAREVEfNHVqsDu7Dtsza5BVY/mIQaVUjPnjvbE8zg83j3OHRMy5AfQDui7ISw9BUbQT8qIMiNXN1kvdI7rDgCXQ+SXxWEEiIhoyDAaIiIisUOsMOFzYgO1XanC0uMniEYMiADeHuGNFnB/mRXvDUc4jBulaIlUz5CU9JwkchEhneSAlAGh9JhnDgIgl0HtE81hBIiKyCQYDREREVxEEARcrW7EjsxYZOXVoU1s+YjDM0wHL4vywNNYX/q5KG3dJI524vQryol1QFO6CrPIYRAbLzyNBJIY2cBo04d3HCroE2bhTIiIiBgNEREQAgPLmLqRn1mJHVg3Km1UWa9yUUiye4Itl8X6I8+MRg3QtSVO+aXhgrycJSBTQjJttHB4YtgiCg6ftmiQiIrKAwQAREdmtVpUWe3LrseNKDS5UtlqskYpFmBXpheVxfrg13INHDNL3BAOkNeehKNoFeeFOSJsLrJYa5K7QhC3oPklgLiDnQEoiIho5GAwQEZFd0ekNOFrchB2ZNThU0ACt3nxuAAAkBLhgebwfFo73gZsDJ8BTN70WsspjxnkBRbsg6aixXuroB03EYqgjlkAbeAsgkduwUSIior5jMEBERGOeIAjIrG7Djsxa7M6pQ3OX1mJdoKui+4hBP4R48Gx46qbpgLx0v3FlQMk+iNUtVku/P0lgMXR+yTxJgIiIRgUGA0RENGZVtqiwM6sWOzJrUNJkeRK8s0KCheN9sCzOD4lBrhBzbgABEHXWQ1GcAXnRLsjLvoNIr7Zaq/VNhDpiKTThi6H3jLZhl0RERIODwQAREY0pbSod9ubWYUdWLc6VW35nVyICpod7YlmcH2ZFeEIp4xGDBIhbSrrnBeyCrPoURILBYp0glkIbOB3qiMXQhKfA4Bxo406JiIgGF4MBIiIa9bR6A44VNyG9e26AxsrcgFg/ZyyL80PKBB94OnK/t90TBEjrrxhPEijaCWlDtvVSqQM0ofOgDl8MTegCCEp32/VJREQ0xBgMEBHRqCQIAq70zA3IrkWLyvI58X4uCiyN9cWyOD+EeznauEsacQw6yCpPQF60C4rCXZC0V1gvVXpCHb4ImvAl0IybCUg5d4KIiMYmBgNERDSqVLR0IT2zFulZtSi1MjfASS7BgvHeWBbnh+RgN84NsHfaTsjLDkJRuAvy4j0Qq5utlupdxhmPFIxYDK3/FEDMbSZERDT2MRggIqIRr1WlxZ7ceqRn1uB8RavFmp65AUtjfTE70otzA+ycqKsBiqIMyIt2Q152sPfhgd4TjccKhi+G3isWYJBERER2hsEAERGNSBqdAUeKGpGeVYvDhQ3Qcm4AXYe4pRiKot3XHx4okkAbOA2acGMYYHANtnGnREREIwuDASIiGjEEQcCFilakZ9ViT24dWq3MDQhwNc4NWBrrhzDODbBfggBp3SXjvICiXdcZHqiEZtwc47GCYQsgKD1s2CgREdHIxmCAiIiGXXFjJ9KzarEzswaVrZaXfDsrJFgw3gfL4nyRFMS5AXZLr4Ws8jgURTshL9oNSXuV1VKD0gPqsBRoIhZDEzwLkHF4IBERkSUMBoiIaFg0dGiQkVOH9KxaZFa3WayRikWYEe6JZXG+mBHhBYVUbOMuaSQQadogLzkAedEuyEv2QayxPGcCAPSuocYjBSNSoPWfDIj5pw4REdH18F9LIiKymS6tHgfzG5CeVYMTxU2wMjYAkwJdsSzOFwvG+8DdQWbbJmlEEHfUQF6UYQwDyo9AZNBYrdX6JEATngJ1xBLoPSdweCAREVE/MRggIqIhpTMIOF3ahPSsWuzPq0eX1vJAuBAPByyN9cWSWF8Eu3PJt90RBEia8iAv2g1F4U7Ias9bLxVLoQ2cDnV4CjThKTC4BNmuTyIiojGIwQAREQ06QRCQXduOnVm12JVdh4YOy+/2ejjIkDLBB0vj/BDn5wwR3+m1LwY9pNVnoCjaBXnRLkhbiq2XypygCZlnnBcQMg+C0t1mbRIREY11DAaIiGjQVLR0YVdWHdKzalDc2GWxRiEVY26UF5bG+mFaqDukEs4NsCvaLsjLDhlXBpTsgbirwWqp3tEXmrBF0ISnQDNuJiBR2LBRIiIi+8FggIiIbkhzlxZ7cuqwM6sWFyotD4UTi4ApIe5YGuuHudFecJLznx97IupqgLx4DxRFuyEvOwiRTmW1VucRbZwXEL4YOr8kQMTgiIiIaKjxLzMiIuo3lVaP7wobsTOrFkeLGqEzWJ4iOMHXGUtifZEywQc+zny3156Im4ugKNoNRdEuSKtPQyRYni0hQARdwGTTsYJ69wgbd0pEREQMBoiIqE/0BgGny5qRnlWLA3n16NDoLdYFuCqwpHuIYISXk427pGEjGCCtOWdcFVCUAWlTrvVSiQKacbONKwPCFkJw9LFho0RERPRDDAaIiMgqQRCQU9uO9Kxa7M6uQ72VIYKuSikWjvfB0lhfTApyhZhDBO2DTgV5+RHjkYLFeyDprLVaalB6QBO2yHiSwLjZgMzRho0SERFRbxgMEBGRmfLmLuzKrsXOrFqrQwTlEhFmRXphyQRf3BruCbmUe8HtgUjVBHnxXuNJAqUHIdJ1Wq3Vu4ZCHb4YmogUaP0nA2L+2UFERDQS8V9oIiICADR1apCRU4+dWbW4VGV5iKAIwOQQdyyJ9cX8aG84K/jPiD0QtxR3bxHYDVnVKYgEy9tIAEDrmwRN+GKow1Og9xwPcPUIERHRiMe/6IiI7FiXVo+D+Q3YmVWL4yVN0FsZIhjTM0Qwxge+LhwiOOYJBkhrzkNenAFF0W5IG3Osl4rl0ATPgCZ8MTThC2Fw8rdho0RERDQYGAwQEdkZnd6AE6XN2JlVi4P59ejSWp4WH+iqwGIOEbQfpnkBu7vnBdRYLTUo3KAJXQB1eAq0IXMhyJ1t2CgRERENNgYDRER2QBAEXKpqw86sWuzJqUNTl9ZinZtSioUx3UMEA10h4jLwMU3U1Qh58R4oind3zwuwPE8CAPSuIcbBgeEp0AZM5bwAIiKiMYT/qhMRjWFFDZ3YmVWDndl1qGxRWaxRSMWYE+mFJbG+mB7mAamEQwTHMklzIeRFu41bBKpPQyRYXjEC9MwLSOmeFxDDeQFERERjFIMBIqIxprZNbTpRILeuw2KNRARMCfXA0lhfzInygpOc/xyMWQY9pDVnTcMDpc0FVksFicI4LyAshfMCiIiI7Aj/EiQiGgNaurTYl1ePXdm1OFvWAssjBIGEABcsifXFgvE+8HKS27RHsiFtJ+RlhyAvyoCiZA/EXQ1WSw0Kd2jCFhq3CYybA8g5T4KIiMjeMBggIhqlVFo9vitsxK6sWhwpaoTOyokCYZ4OWBLri8UTfBHs7mDjLslWxB01kBfvMQ4PLD8MkV5ttVbnFtZ9isAiaP0nc14AERGRneNfAkREo4jOIOBkSRN2ZdfiQF4DOrWWz5P3dZYjZYIvlkzwxXhfJw4RHIsEAZLGbCiKMiAv2g1Z7XnrpRBB53+TcVVA2CLoPaI5L4CIiIhMGAwQEY1wgiDgYmUrdmXX9XqigItCigXjvbEk1hdJQW6QiPnCb8zRayGrOgl50S4oijIgaSuzWipIldCMmwNN2CKowxZAcPSxYaNEREQ0mjAYICIaofLrOrAzuxYZ2bWobLW8LFwhFWNWhBeWxPpgepgn5FKeKDDWiNQtkJceMG4RKNkPsabVaq3BwQfq8IXG4YHjZgJSbh0hIiKi62MwQEQ0glS0dGF3dh12ZdeioL7TYg1PFBj7xC0lUBRnQF6UAVnVCYgMOqu1Os8Y0xYBnV8SIGI4RERERP3DvyaJiIZZQ4cGe3PrsDOrDpeqrL8bnBDgiiWxPjxRYCwSDJDWnDPOCyjOgLQxx3qpSAJt4C3QhC+COmwRDG6hNmyUiIiIxiIGA0REw6BdrcP+vHrszq7DydImWDlQABFejlgS64uUCT4IcuOy8DFF2wl52XeQF++GongvxF31VksNchdoQuZBE74ImpB5EJTutuuTiIiIxjwGA0RENqLS6nG0qBE7s+twpLABGr3lNCDAVWE6USDKh2fKjyXijmrIi/ZAXpxx3SMF9S7joA5fBE1YCrSBUwEJV4kQERHR0GAwQEQ0hHQGAadKm7Aruw4H8urRobF8vKCHgwwLY3yweIIPJgW68njBsUIQIKnPhKJ4N+TFeyCrvdBrudYvGZqwFKjDF0LvOYFHChIREZFNMBggIhpkBkHAxYpW7Mquxd7ceqvHCzrJJZgb7Y3FE3wwJcQDUh4vODbo1ZCVH4Wi2LgyQNJeabVUkCqhCZ4FTdhCaMIWwuDkZ8NGiYiIiIwYDBARDQJBEJBb24Fd2bXYnVOHmjbLS8TlEhFmRHhhyQQf3BruCaVMYuNOaSiIuhogL94LRXEGZGWHINZ2WK3VO/pBE7YAmvAUaIJmADLOjiAiIqLhxWCAiOgGlDR2mo4XLGnqslgjEQFTQjyQMsEH86K94azgr95RTxAgacqDvDgDiqIMSKvPQAQrEyQB6LziuucFLILOdxKPFCQiIqIRhX+dEhH1U3WrChk5ddidXYfs2nardYmBrkiZ4IuFMd7wdOTguFFPr4Gs8qQxDCjeA0lridVSQSyHNng61GGLoAldCINrsA0bJSIiIuofBgNERH3Q2KnB3tx67M6uxfmKVqt1432csHiCLxZN8EGAq9KGHdJQEKmaIC/ZB3nxHshLD0CsabNaa1B6QBO2EOqwhdCOmwNB7mzDTomIiIgGjsEAEZEVbSodDuTXY3d2HU6VNsHK6YII8XBASowPUib4ItzL0bZN0uASBEiaCyAvyjCeIlB9CiLBYLVc5xFtDAPCU6DzuwkQc2YEERERjT4MBoiIrtKl1eO7ggZk5NThSFEjtFbSAF9nOVIm+GLxBB/E+DrzeMHRTK+FrOqkcVVAcQakLcVWSwWxFNqAadCEL4I6bCEMbmE2a5OIiIhoqDAYICK7p9EZcKy4Cbuza3GooAEqneV3iD0cZFgw3huLJ/hiUpArxAwDRi2Rqhny0v3fbxFQt1itNSjcoAmZB034ImhC5kJQuNmwUyIiIqKhx2CAiOySziDgTGkzdmXXYn9+PdrVeot1TnIJ5kZ7Y/EEH0wJ8YBUzDBgtJI0FZhWBciqTkEkWP6eA4DOPQKa0IXQhC+CNmAKIOY/l0RERDR28S8dIrIbBkHAhYpW7M6uxd7cejR1aS3WKaRizI70QkqMD6aHe0Ih5dFyo9I1WwT2QNpSZLVUEEmgDZgCTdgiaMIXQe8eYcNGiYiIiIYXgwEiGtMEQUBmTTt2Z9diT04dats1FuukYhFuDfdESowPZkV6wVHOIXKjkfEUgau2CGisnyBhkLtCEzoPmrCFxi0CSg8bdkpEREQ0cjAYIKIxRxAE5Nd3ICOnDruz61DRorJYJxYBk8e5Y/EEX8yN9oKrUmbjTumGCQIkTfmmVQHXPUXALdy4KiBsAbQBUwEJv+dEREREDAaIaMwobuxERnYdMnLqUNTYabUuOcgViyb4Yn60N7yc5DbskAaFXgNZ5UnIizOgKN4DSWuJ1VJBJIE2cKppXgC3CBARERGZYzBARKNaRUuXKQzIreuwWhfn74KUGB8sGO8Nf1elDTukwSDqavh+i0DZQYg1bVZrTacImLYIuNuuUSIiIqJRiMEAEY06tW1q7Mk1hgGXq6y/QIzydkLKBB8sivFBsLuDDTukGyYIkDRmQ168F4riDEirz0IEwWq5ziMKmtAFxlME/CfzFAEiIiKifuBfTkQ0KjR0aLA3tx57cmpxvqLV6kvEEA8HpMT4YNEEH0R4Odm0R7pBOhXkFUchL94LefEeSNorrJYKYim0AdOgCV8EdegCGNzDbdgoERER0djCYICIRqzmLi3259Vjd04dzpY1w2AlDQh0VWBhjC9SJvhgvI8TRCKRbRulARN3VENess8YBpQdgkjXZbXWoPTo3iKwCJqQORAUrjbslIiIiGjsYjBARCNKm0qHA/n1yMipw8nSZuitpAE+znIsijFuE4j3d2EYMFoIBkjrLnWfIrAXsrqLvZbrPGOgCVsIddhC6PxuAsQ8RpKIiIhosDEYIKJh16HR4VBBAzKy63C8pAlaveUwwNNRhgXjjWFAYpArxAwDRgdNB+Tl3xnDgJJ9kHTWWi0VxHJog6dDHboQmrCFMLiOs2GjRERERPaJwQARDYsurR7fFTRgT249jhY1Qq2zfPa8m1KKedHeWBTjg5vGuUMqZhgwGohbS42DA0v2QFZ+DCKDxmqt3tEXmtD50IQtgCZ4NiDnbAgiIiIiW2IwQEQ2o9LqcbS4CRnZdThc2ACVlTDASS7B3O4wYFqIO6QSsY07pX4z6CCrPt09OHAvpE25vZZrfSYZg4CwhdD5JAAifo+JiIiIhguDASIaUhqdAceKm7Antw6H8hvQqdVbrHOUSTAr0hOLYnwxPcwDcilfKI50IlUT5KUHjFsESg9ArG6xWitIHaAZN9sYBoTOh8HJ34adEhEREVFvGAwQ0aDT6g04WdKMjJxaHMhvQIfGchigkIoxK8ILi2K8cWu4J5QyDpYb0QQBksYcyEv2QlG8F9Lq0xAJlld9AIDeZRw0YQugDl0AbdB0QKq0YbNERERE1FcMBohoUOj0Bpwqa8aenDocyG9Aq0pnsU4uEeHWcE8sivHBrEgvODAMGNl0Ksgrjhq3CJTshaSt3GqpIBJD5z8Z6rAF0IQuhN5zPMABkUREREQjHoMBIhownUHAme4wYH9ePVqshAEyiQjTwzyxMMYbsyK84Kzgr56RTNxeCXnxPshL9kJefhgiXZfVWoPCDZqQudCELYQmZC4EpYcNOyUiIiKiwcC/zomoX/QGAWfLm7Enpx778urR3KW1WCcRi3BLqAcWxfhgdqQXXJT8dTNiGfSQ1pz7fotAQ2av5TrPmO5TBBZC638zIOb3loiIiGg0419zRHRdeoOA8xUtyOheGdDYaSUMEAFTQj2waLwP5kR5wc1BZuNOqa9EqmbIyw4atwiU7odY1WS1VpAooA2aDnXYQmhCF8DgOs6GnRIRERHRUGMwQEQW6Q0CLlS2mFYGNHRYPodeLAImj3PHwhgfzIvyhrsjw4ARSRAgacw1bg8o3gtZ9WmIBMtDIQFA7+QHTegC4/+NmwXIHG3YLBERERHZEoMBIjIxCAIuVLRiT04d9uXVo76XMOCmce5YNN4bc6O94ekot3Gn1Ce6LsjLj0Jesu/6gwMhgs4vufsUgYXQe8dxcCARERGRnWAwQGTnDIKAixWt2JNrDAPq2i2HASIAycFuWBjjg/nR3vByYhgwEonbKkxBgHFwoMpqrUHuCk3IHGjCFkATMg+Cg5cNOyUiIiKikYLBAJEdMggCLlW2mmYG1PYSBiQFu2HheG/Mj/aGt7PCto3S9Rl0kFWfMW0RkDbm9Fqu8xgPTdh8aEIXQOs/GZBw6wcRERGRvWMwQGQn+hUGBLliwXgfzB/vDR+GASOOqKsR8tL9xpUBpQcgVrdYrTUNDgxdAE3YAhhcQ2zYKRERERGNBgwGiMawnjBgT2499uXWWQ0DgKvCgGhv+LowDBhRBAHS+ivfrwqoOQcRBKvleucAaEKMxwlqgmdwcCARERER9YrBANEYc/XMgN5WBgBAYqArFnTPDPBjGDCiiDTtkJUd6p4XsB+SzhqrtYJIDJ3/zcZVAaHzofeK5eBAIiIiIuozBgNEY0DPaQJ7rzNAEAAmBbqaBggyDBhBBAGS5kLjqoCSfZBVnoDIoLVablC4QxM6z3icYMgcCEoPGzZLRERERGMJgwGiUUpvEHChsgV7c+p7PVoQMIYBC8Z7Y8F4H4YBI4lOBVnlcciL90JRsg+S1pJey7Xe8cYgIHQ+dH7JgFhio0aJiIiIaCxjMEA0iugNAs5XtGBvrjEMaOhDGDA/2hv+rkobdkm9+f44wX3dxwl2Wa0VpI7QjJtlOk7Q4Bxgw06JiIiIyF4wGCAa4XQGAefKm7E3tx778+rR2Gl5ebkIQCIHCI48ei1kNWeMQUBfjhN0CzcGAaHzoQ2cBkj4fSQiIiKiocVggGgE0hkEnClrxt7cOhzIa0BTl/UwgEcLjjyijlrISw8Yw4CyQxBrWq3WCmI5tEHToQmdD03oPOjdI2zYKRERERERgwGiEUOnN+BUmXFlwIG8erSodBbrRACSg92wYLw35kUzDBgRDHpIay+YtgjI6i72Wq53DjDNCtAEz+RxgkREREQ0rBgMEA0jrd6AkyXGlQEHCxrQaiUMEIuAm4LdMH+8D+ZFe8PbSW7jTumHRKomyEsPGsOA0gMQqxqt1goiCbQBk6EJmQdN2ALoPSfwOEEiIiIiGjEYDBDZmFpnwImSJuzrDgPa1XqLdWIRcNM4dywc7425Ud7wYhgwvAQB0vorplUB0pqzEAkGq+UGB2/jcYIh86EJmQ1B4WbDZomIiIiI+o7BAJENqLR6HCtuwt7cOhwubESHxnIYIBEBU0I8MH+8N+ZGecHDkWHAcBKpWyErOwR56X7ISw5A0lljtVaACDrfROP2gLAF0PkkACKxDbslIiIiIhoYBgNEQ6RTo8fRokbsza3HkaIGdGktv7ssEYswLdQdC6J9MDvKC+4OMht3SiaCAEljzvezAqpPQ2SwvL0DAAwKN2hC5hpXBoybC8HR24bNEhERERENDgYDRIOoXa3D4cJG7M2tw7HiJqh1lsMAmUSEaaEeWDDeG7MjveCqZBgwbDQdkJcf7p4VsB+S9spey7XeE41BQOh86PySATF/jRIRERHR6Ma/aIluUKtKi0MFDdibW48TJU3Q6gWLdQqpGNPDjNsEZkV4wVnBH79hIQiQNBd0rwrYD1nlCYgMGqvlBpkztONmdZ8iMBcGJ38bNktERERENPT4yoRoAJo6NTiQ34B9efU4VdoMvcFyGKCUijEzwhPzx/tgRrgnHOUSG3dKAABtJ+QVRyEv2Q95yT5I2sp6Ldd5xphWBWj9JwMSznogIiIiorGLwQBRH9W3q7E/vwH7cutwtrwFVrIAOMokmBVpDANuDfOAUsYwwOYEAZKWoqtWBRyHSK+2Xi51gCZ4lnFwYOg8GFyCbNgsEREREdHwYjBA1IvqVhX25dVjX249Lla2wkoWABeFFLO7w4BpoR5QSDmN3ua0XcZVAaX7IS/ZD0lrSa/lOvfI7uME50EbOA2QKm3UKBERERHRyMJggOgHypq6jGFAXj0yq9us1rk7yDAnygvzo70xJcQdMgnDAFuTNBcatweU7oes4th1VgUooQm61bgqIGQeDG6hNuyUiIiIiGjkYjBAdk8QBBQ1dmJfrjEMyKvrsFrr5STHvCgvzB/vjeRgd0jFIht2StB1QV5xDLLuMEDaUtx7uVuYKQjQBt0CSB1s0ycRERER0SjCYIDskiAIyKltN20TKGnqslrr56LA/GhvzI/2RkKgKyQMA2yqX6sCJIruVQHGLQIG93AbdkpERERENDoxGCC7YRAEXKpsxb68ehzIq0dlq/UXmMHuSlMYEOfvApGIYYDN9HNWgN411HSCgCZoOlcFEBERERH1E4MBGtN0BgHnypuxL7ceB/IbUN9h/bz6cC9HUxgQ7ePEMMBWBMG4KqC0Z1XAdU4QkCigDboFmhBjGKB3j7Bhs0REREREYw+DARpzNDoDTpU2Y19eHQ7mN6BFpbNaO8HXGfO6w4AwL0cbdmnntJ3GVQHdWwQkraW9lptWBYTMgyboVkDGVQFERERERIOFwQCNCV1aPY4VNWJfXj0OFzaiQ6O3WpsQ4Ir5470xL9oLQW58gWkTggBJUz7kpQe+XxVgsL56w7gqYHr3qoB5XBVARERERDSEGAzQqNWm0uG7wgbsz6vHseImqHUGi3ViEZAc7Ib50d6YG+UNXxeFjTu1TyJNO2TlR7q3CByApK2813qdWzg0IXOhDZ3HWQFERERERDbEYIBGlYYODQ4WGMOAU6XN0BsEi3VSsQjTQj0wL9oLsyO94OEot3GndkgQIGnMhryke1VA1SmIDFrr5VKl8QSB7lUBBrcw2/VKREREREQmDAZoxKtuVZlOEjhf0QrLUQCglIpxa7gn5kV7Y2aEJ5wVfHoPNZG6FbLy77q3CByApL2q13qdeyQ0IXOhCZ0HbeAtgFRpo06JiIiIiMgavnKiEam4oRP78+uxP68eWTXtVuucFRLMjPDC/GhvTA/zgFImsWGXdkgQIK2/AlnpAchL9kNWfRoiwfo8B0HqAE3wzO7BgXNhcA2xYbNERERERNQXDAZoRBAEAdm17difV48DeQ0oauy0WuvpKMOcKC/Mi/bG5HHukEnENuzU/ohUTZCXfWfcHlB6EJLO2l7rdR7jTScIaAOnABLOdCAiIiIiGskYDNCw0RsEXKhswf68BhzIq0d1m/Wz6/1dFJgX7Y150d6YFOgKiVhkw07tjEEPae0F0/YAae15iATLgx0BwCBzhnbcTOMWgZB5MLgE2bBZIiIiIiK6UQwGyKY0OgNOlTZjf349DuU3oKnL+nC6UA8HUxgQ6+cMkYhhwFARddRCXnbIeIJA2SGIVU291uu84kzbA7T+kwGJzEadEhERERHRYGMwQEOuQ6PD0aImHMirx5GiRnRorO9Jj/VzxtwoYxgQ7uVowy7tjF4LWc0ZyEsOQFZ6ALL6y72WGxRu0ATPMg4NDJkDg5O/jRolIiIiIqKhxmCAhkRzp/b/t3fn4TGe6x/Av5PJTDaCLFQkiGUGycGxBCHLiQiHQyjiWIP2qBY9rdovVC1HVaunOLSK1tbSSGlRtQdtlNZSRZOIBBMqQmRfZibz/P6ImZ+RzGQpk2W+n+vqH3nfe973mWduaeZ+nwWnbjzEicQHOHfrEdRFpe8lIAHQ0bMegh+vGdDYmavUPy82WSmQq4qnB8hUP8BGY3pRRwEJtA07GHYQ0DbsANjw1wURERERUW3Ev/TpmfkjqwAxiQ8fbyuYCZ2JfQVtbSTwa1Yfwa3cENjSFa5Ocss21Fpo8yG7exby2yeL1wp4dN1suM7BDeqmQcXFAK8gCAcXCzWUiIiIiIiqEgsDVGlCCCQ9zENMYvFOAnH3TT+BdpDZoKe3C4JbuaFnCxfUsWPqPXNCQJpxw7BooOzOGUiKTC/oKCRSaF7oAk3TYKibBkHr7gtIuMMDEREREZG14bczqhCdELj6R3ZxMSDxIW4/yjcZW8/eFoEti6cI+DVrADtbful81iTqbMhSfjCMCpBmp5iNL6rj8Xj3gGBoPHtB2DlbqKVERERERFRdsTBAZdIU6XBelYGYxIc4mfgQD3LVJmMb1bUzrBfQoUk92HJbwWdL6GD74GrxooGqGMjunYdEpzUdbiOHxqObYQeBogatAe7uQERERERET2BhwIzLly/j7bffRnJyMpRKJZYvX44WLVpUdbMsIk9dhNjkdMQkFu8kkFNoeicBb1dH/K2VK4Jbu6FNQ24r+KxJ8h5ArjpZPCpAdQo2+Q/MxmvrtygeEeAVBHWTHoCMuzsQEREREZFpEiGEiSXirFthYSF69+6NGTNmoH///tiwYQNiYmKwe/fuSl0vLS37Gbfw2csp1OJYQhpiEh+a3UkAAHwb10VwKzcEt3JFMxd+8XymHm8lKHs8PUCW9pvZcJ3MCRrPXo+nCARB59zUQg0lIiIiIqLqzN29brniOGLAhJ9++gl16tTB4MGDAQCvvvoqtmzZgvj4eCiVyqpt3HNQoCnCqK3n8UdW6YvVSW0k6OJVD8Gt3BDUyhXudews3MLazSbrtmGdAFnKj2a3EgQAjZuPYdFAzQtdACl3diAiIiIiosphYcCE5ORko2kDUqkUXl5euHHjRq0sDKgy8ksUBRxkNvD3dkFQK1f08nZFXXumyzOjyYP8zhnIHu8gYJuZbDZcZ+8CtVfg460EAyGcGlqooUREREREVNvxm54JeXl5cHBwMDrm4OCAgoKCKmrR89XSzQlDOzTGxZRMwzSBrk3rw14mreqm1Q5CQPrw98frBJyE7O45SHSmF3EUEim0L3SGumkQ1E2DoXX/C7cSJCIiIiKi54KFARNKKwLk5+fD0bF2zqe3kUgwJ7R1VTejVpHkp0OuOlVcCLh9CtK8VLPxRXU9ofYKgrpZMDRNenIrQSIiIiIisggWBkxo0aKF0UKDRUVFuH37ttXsSkCVUKSBLPVC8aKBqpOwvX8ZEphewFHY2kPt0ePxWgHBKKrfglsJEhERERGRxbEwYEK3bt2QkZGB6OhoDBw4EBs2bICHhwdat+ZTdfp/Npm3DFsJyu7EwkZtfvcJrYvy8e4BwdA07grY2luopURERERERKWr1dsVzp8/H7du3cK2bdtKnFOpVFixYgXOnTsHAAgODsacOXPg4uJiiLly5QrefvttJCUloU2bNli+fDmaN29eqbbUhO0KqRzUuZDfiYVcFQPZ7ZOwzbxpNlxnV//xooFB0HgFQlensUWaSUREREREZPXbFUZFRSEqKgp+fn4lzj169AiRkZFQq9V4+eWXUVRUhE2bNiE+Ph5RUVGQy4u3fvP19UV0dLSlm07VidDB9sHV4ukBt2Mgu3ceEp3GdLhECu0LnYrXCmgaBK17e8CGCzgSEREREVH1VesKA0VFRVi/fj3Wrl1rMubzzz/HvXv3sG/fPrRs2RIA0KFDB0yYMAF79+5FRESEpZpL1ZBNbipkqlPFOwiknIZN/kOz8UV1vR7vHhAETRN/CLt6FmopERERERHRn1erCgOFhYUYPnw44uPjMXjwYJw5c6bUuAMHDsDPz89QFAAAf39/eHt748CBAywMWBttAWR//Az57RjIVadg+/B3s+HC1gHqJv7FhYCmwSiq581FA4mIiIiIqMaqdYWBnJwcfPjhh+jfvz9CQkJKxGRmZkKlUqFv374lzvn4+CAmJsYCLaUqJQSkj64XjwhQnYTs7k+QaAvMvkTj5gtN00CovYKgadwFkNpZqLFERERERETPV60qDNSpUweHDx+Gra3pt5WaWryXfKNGjUqcc3d3R05ODrKzs1G3bvkWaaCaQVLwCHLVacNWgtLce2bjdQ5uhkUD1V6BEI7uFmopERERERGRZdWqwoCNjQ1sbGzMxuTm5gIAHBwcSpyzsyt+CpyXl8fCQE1XpIYs9QJkt09BrjoJ2/uXIYHpDTiEjRyaxl0eFwKCUeTWFpCYzyUiIiIiIqLaoFYVBspDp9OVGVNWcYGqISEgzUx+PCLgFGR3YmGjyTX7Em2D1lB7BULjFQh1kx6AzNFCjSUiIiIiIqo+rK4w4OTkBKB4PYKn6Y/pY6h6kxRkQHbnR8hvn4JcdQrSbJXZeJ1dPag9A6B5PD1AV7eJhVpKRERERERUfVldYcDDwwMAkJaWVuLc/fv34ezsDEdHPjmuloo0sE29CLnqJOSq07C9fwkSYXoEiLCxhaZR5+JFAz0DoG3YAbCRWrDBRERERERE1Z/VFQacnZ3h6emJq1evljh37do1+Pr6VkGrqFRCwCbzJuSq4hEBspQfYaPJMfsSbb3m0HgVjwjQePpDyLlWBBERERERkTlWVxgAgLCwMGzduhU3btxAy5YtAQCxsbFITk7GSy+9VMWts24Vnh4gd4bG0x9qr2CovQKgq9fMQi0lIiIiIiKqHayyMPCvf/0L33zzDcaPH4+JEyeisLAQGzduRLt27RAeHl7VzbMuRZri3QMejwqwvf+r+ekBEim0jf5q2EqweHqAVaYxERERERHRM2GV36hcXFywfft2LF++HKtXr4a9vT169+6NmTNnQi6XV3Xzajf97gGqU5Df1u8eYH56QJFzs8eFgEBomvhD2NWzUGOJiIiIiIhqP4kQwvTm7vTMpKVlV3UTqoyk4BFkKT8aFg2UZqeYjdfJ60Lj2bO4GOAZAF19bwu1lIiIiIiIqPZwdy/fmmtWOWKAnrMiNWT3zj8xPeAyJDBdfzKaHuAVCG2jjpweQEREREREZCH89kV/nhCQPrpevHOA6hTkd36CRJtn9iWcHkBERERERFQ9sDBAlSLJewB5ymlDMUCam2o2XmdXr3h6gGcgdw8gIiIiIiKqRlgYoPLR5kN29xzkqlOQq07D9uE1s+HCxhaaRp2h8Qoonh7QsANgI7VQY4mIiIiIiKi8WBig0gkdbB9cfbxOwGnI/vgZkqJCsy/RNmgFtWcANF6B0DTpASGvY6HGEhERERERUWWxMEAGksJM2N04CFnKachVp2FTkG42XmfvArVXADSexaMCdHU9LNRSIiIiIiIielZYGKBimjw0+LI3pLn3TIYIqR00jbsWFwO8AqF18wEkNhZsJBERERERET1rLAwQAECaebPUooDWtR3Uj9cJ0DT2A2QOVdA6IiIiIiIiel5YGCAAQJFrG+R1eBmyO2dQ5OZTXAzw7AXh6F7VTSMiIiIiIqLnSCKEEFXdCGuQlpZd1U0gIiIiIiIiK+LuXrdccZwgTkRERERERGTFWBggIiIiIiIismIsDBARERERERFZMRYGiIiIiIiIiKwYCwNEREREREREVoyFASIiIiIiIiIrxsIAERERERERkRVjYYCIiIiIiIjIirEwQERERERERGTFWBggIiIiIiIismIsDBARERERERFZMRYGiIiIiIiIiKwYCwNEREREREREVoyFASIiIiIiIiIrxsIAERERERERkRWTCCFEVTeCiIiIiIiIiKoGRwwQERERERERWTEWBoiIiIiIiIisGAsDRERERERERFaMhQEiIiIiIiIiK8bCABEREREREZEVY2GAiIiIiIiIyIqxMEBERERERERkxVgYICIiIiIiIrJiLAwQERERERERWTEWBoiIiIiIiIisGAsDRERERERERFaMhQEiIiIiIiIiK8bCAJGVOH36NEaNGoUOHTrgr3/9K8aPH49Lly4ZxahUKkydOhV+fn7w8/PDrFmzkJ6eXjUNJiqHuLg4+Pr6Ys2aNUbHmctUk6Snp2P+/Pnw9/dHp06dMGbMGFy4cMEohjlNNcWVK1cwYcIEdOzYEZ06dcLkyZORlJRkFMN8pupu/vz5GDt2bInj5c3dmpjjEiGEqOpGENHzdfbsWURGRqJ169YYOnQotFotvvjiC9y/fx87duxAhw4d8OjRIwwdOhRqtRrjxo1DUVERNm3ahCZNmiAqKgpyubyq3waREa1Wi+HDh+PatWuYOnUqpk2bBgDMZapRcnJyMHz4cNy/fx/jx4+Hs7MzduzYgXv37iEqKgpKpZI5TTVGUlIShg4dCgcHB4wfPx4A8Nlnn0EIgW+++QaNGjViPlO1FxUVhfnz58PPzw/btm0zHC9v7tbYHBdEVOsNHDhQBAcHi7y8PMOxtLQ00bVrVxEZGSmEEGLVqlWibdu2IjEx0RDz448/CoVCIXbt2mXpJhOVae3atcLHx0coFAqxevVqw3HmMtUkq1atEkqlUpw7d85w7P79+6J9+/ZixowZhhjmNNUECxcuFAqFQly9etVw7NdffxUKhUK8++67QgjmM1VfWq1WrFmzRiiVSqFQKMSYMWOMzpc3d2tqjnMqAVEtl5mZiYSEBPTr1w8ODg6G425ubujatathOsGBAwfg5+eHli1bGmL8/f3h7e2NAwcOWLrZRGbFx8dj/fr1eO2110qcYy5TTSGEwJ49exAcHIyuXbsajru7u2PWrFmGY8xpqilSUlLQoEEDtGvXznCsffv2qF+/PhISEgAwn6l6KiwsxJAhQ7BmzRqEh4ejUaNGJWLKm7s1NcdZGCCq5erUqYPvv//eMKTvSY8ePYJUKkVmZiZUKhV8fHxKxPj4+ODKlSsWaClR+Wi1WsydOxf+/v4YNGiQ0TnmMtUkKSkpSE1Nhb+/P4DiQkFubi4AYPTo0YiIiGBOU43SrFkzZGZmGs2lzsjIQHZ2Ntzd3ZnPVG0VFhYiJycHH374IVasWAFbW1uj8+XN3Zqc4ywMENVyUqkUzZs3L1H5jIuLw4ULF9CpUyekpqYCQKnVUXd3d+Tk5CA7O9si7SUqy6effopbt25h8eLFJc4xl6kmuXXrFoDiEVwrV65E165d0alTJ/Tp0wfHjx8HwJymmuXll19G48aNMX36dMTFxSE+Ph5vvfUWbG1tMWbMGOYzVVt16tTB4cOH0b9//1LPlzd3a3KOszBAZIVyc3Mxe/ZsAMArr7xieEL15FQDPTs7OwBAXl6e5RpIZML169fxv//9D7Nnz8YLL7xQ4jxzmWqSrKwsAMBHH32EEydOYN68eVixYgXs7e0xZcoUxMbGMqepRvHw8MCkSZPw888/Izw8HIMGDcKZM2fw/vvvw9fXl/lM1ZaNjU2JUQJPKm/u1uQcN/3uiahWys/Px+TJkxEXF4dXX30VXbp0wfnz58t8nY0N64hUtYqKijB37lx07twZERERpcbodLoyr8NcpupCrVYDKB56eujQIdSrVw8AEBISgj59+uCDDz7AvHnzyrwOc5qqi48++gjr1q2Dn58fIiIiUFRUhC+++ALTp0/Hf//7XzRo0KDMazCfqToq798XNfnvEBYGiKxIZmYmXnnlFVy8eBHDhg3DG2+8AQBwcnICUDy/6mn6Y/oYoqqyadMmxMXF4YsvvjDMX9U/cc3Pz0d6ejpzmWoUR0dHAEBYWJihKAAAzs7OCAkJwZ49e5jTVGNkZWVh48aN8PHxweeffw6pVAoAGDBgAIYOHYqFCxdi8+bNAJjPVPOU93dxTf6dXT3LFUT0zD18+BBjx47FxYsXMWLECCxdutRwzsPDAwCQlpZW4nX379+Hs7Oz4Q9Yoqpy+vRpaDQaDB8+HD169ECPHj0wZMgQAMVFgx49eiAlJQUAc5lqBv0cVBcXlxLnXFxcIISAq6srAOY0VX83b96EWq3GP/7xD0NRAABkMhkGDRqEhw8fQggBgPlMNU95/1auyX9Tc8QAkRXIycnBxIkTER8fj/Hjx2Pu3LlG552dneHp6YmrV6+WeO21a9fg6+trqaYSmTR79mzDCAG9Bw8eYObMmQgPD8fgwYPh4+PDXKYao3Xr1pDL5UhMTCxxLiUlBXZ2dnBxcWFOU40gl8sBwPDl/0n64dU6nY75TDVSef9Wrsl/U3PEAJEVeOeddxAXF4dx48aVKArohYWF4cyZM7hx44bhWGxsLJKTk02u0EpkSb6+vvD39zf6r1OnTgAALy8v+Pv7o169esxlqjEcHR0REhKCmJgYXL9+3XBcpVLh+PHjCA4OhlQqZU5TjdC6dWs0bNgQe/bsMRpGrVar8c0336BBgwZQKBTMZ6qxypu7NTXHJaK0sh4R1RoJCQkYOHAg6tati3nz5hkN79MLDw9Henq6YfjfxIkTUVhYiI0bN8LLywu7du0yPAkgqk5SUlLQu3dvTJ06FdOmTQMA5jLVKCkpKRg+fDgAYNy4cZDJZNi6dSvy8vKwe/duNG/enDlNNcaRI0fw+uuvo1WrVhg2bBh0Oh2+/vprXL9+He+99x4GDRrEfKYaISQkBE2aNMG2bdsMx8qbuzU1x1kYIKrlduzYUep+70+Kj48HACQlJWH58uX45ZdfYG9vj8DAQMycORNubm6WaCpRhZVWGACYy1SzqFQqrFy5ErGxsRBCoHPnzpg5cyZat25tiGFOU01x5swZrFu3Dr/99hsAoF27dpg8eTICAwMNMcxnqu5KKwwA5c/dmpjjLAwQERERERERWTGuMUBERERERERkxVgYICIiIiIiIrJiLAwQERERERERWTEWBoiIiIiIiIisGAsDRERERERERFaMhQEiIiIiIiIiK8bCABEREREREZEVY2GAiIiIiIiIyIqxMEBERERERERkxVgYICKiKpWQkID58+cjNDQU7du3h7+/P0aOHIkvv/wS+fn5z+QeJ06cgFKpxJo1a57J9Z61Xbt2oXfv3vD19UVAQAAKCwtLxEycOBFKpRInT540e63c3Fx07NgRAQEB0Ol0z6vJJlX3vra059Efa9asgVKphFKpxIYNG8zGzpo1yxCbm5tb6XtevnwZsbGx5YqdM2cOlEolEhISKn0/IiKyLBYGiIioynzyyScIDw/Hnj170KJFC4wbNw6hoaFIT0/HokWLEB4ejsTExKpu5nN148YNvP3228jOzsaYMWMQEREBOzu7EnHh4eEAgO+++87s9Y4ePYr8/HwMGjQINjb833xtd+TIEZPnNBoNTpw48afvERMTgxEjRiApKalc8aGhoZg6dSpcXV3/9L0BYO/evVAqlXj33XfLjJ0yZQqUSiW2bNnyTO5NRGQtbKu6AUREZJ22bt2KVatWQaFQYPXq1fD29jac0+l0+PLLL7Fs2TKMHj0a+/btQ8OGDauwtc/P77//DiEExo0bh6lTp5qMCwsLw6JFi3D06FGo1WrI5fJS4/bv3w8AGDJkyHNpL1Uf7u7u+O2335CamopGjRqVOB8bG4usrCw4OjoiLy+v0vdJT0+v0OiT0NBQhIaGVvp+T2vTpg0AlDkCITY2FkePHoVCocDo0aOf2f2JiKwBHyUQEZHFpaamYuXKlXBxccHmzZuNigIAYGNjg9GjR2P69OnIyMjA4sWLq6ilz59arQYA1K9f32ycg4MD+vbti5ycHJw6darUmPT0dMTGxsLX1xetWrV61k2laiYkJARCCJOjBg4dOoSGDRvCx8fHwi17tlq2bAmZTGa2MKDVarFs2TIAwMKFC2Fry2dfREQVwcIAERFZ3K5du6BWqzFq1Ci4u7ubjIuMjIS7uzuOHTuG+/fvAwDGjh2LPn364NixYwgODkaHDh0wZ84cw2uuXr2KSZMmoUuXLujWrRsWL15scq2CrKwsrFixwjC/PygoCEuXLsWjR4+M4sq6Z2nu3buH+fPno1evXvD19UVISAjeffddZGZmGmJCQkIwd+5cAMCSJUugVCrx9ddfm7ymfjrBwYMHSz1/8OBBaLXaEqMF4uLi8NZbbyEwMBC+vr7o3Lkzxo4dW2K9Av37vHv3Lt588034+fmhY8eOGDduHC5evFjifpbua/05lUqFqVOnonPnzujatStmzJiBrKwspKenY/bs2fDz80O3bt0wbdo0pKamGl2/PJ9LZfriefWHOT179oSTk1OphYGioiIcO3YMffv2hUQiKfX15cmLOXPmlMjRlJQUs5/Tk2sM/PTTT2jTpg1CQkKM+iIvL8/w3q9du2b2fcpkMrRo0QJpaWlIT08vNWbHjh1ITEzEwIED0bVrV/MdR0REJbCcSkREFnfmzBkAQEBAgNk4mUyGwMBAREdH4+TJkxg+fDiA4ifjM2bMQJ8+fWBvb294Inr58mVERkZCq9WiX79+cHBwwMGDB3H48OES187KysI///lP3LhxA7169UK/fv1w8+ZNbN++HadPn8bOnTvRoEEDQ7ype5bm5s2bGDlyJNLT0xEQEIBWrVrh8uXL+Oyzz3DixAnDtceNG4dz587h2LFjCAwMRPv27dG2bVuT1+3evTs8PDxw/PhxFBQUwN7e3uj8vn37IJPJMGDAAMOxixcvIjIyEvb29ggLC0P9+vWRnJyM48eP4+eff8b27dvRpUsXQ3xOTg5GjRqFunXrYsiQIUhNTcX333+P8ePH48iRI4YpHVXV1zk5ORg5ciQ8PT0xYsQInD17Fvv27UNWVhZSUlIgl8sxdOhQXLt2DYcPH8ajR4+wffv2Cn0uFe2L590fpsjlcgQHB+PQoUPIyMgwGnVy9uxZZGRkoG/fvoiPjy/x2vLmRWhoKLKysoxy1NnZuczPSa979+4YMWIEdu7cifXr12P69OkAgJUrVyIlJQVvvPEG2rVrV+Z7VSqViI+PR0JCArp37250Lj09HWvXroWTkxNmzZpV5rWIiKgUgoiIyMK6desmFAqFyMjIKDN2w4YNQqFQiFWrVgkhhBgzZoxQKBRi5cqVJWIjIiJEu3btxLlz5wzH7ty5IwICAoRCoRCrV682HF+4cKFQKBQiKirK6Br79+8XCoVCzJ0713DM3D1Lo4/fu3ev0fHVq1cLhUIh5syZYzgWHR0tFAqF2LZtW7mu/cEHHwiFQiG+//57o+MqlUoolUoxZcoUo+MTJkwQPj4+Iikpyej4zp07hUKhEIsXLy7R7mnTpgmNRmM4vmrVKqFQKMSGDRsMx6qir/XnXn/9daHT6YQQQhQUFIju3bsLhUIhJk6cKLRarRBCCJ1OJ4YOHSoUCoW4d++e0evL87lUpC+eZ3+URt/e48ePi++++04oFAoRHR1tFLNw4ULRs2dPUVRUZHgvOTk5hvMVyYvSctTc5zR79myhUChEfHy8EEKI7OxsERwcLHx8fERiYqL46aefhFKpFBEREYbPqywbN24UCoVCbNmypcS5+fPnC4VCITZv3lyuaxERUUmcSkBERBaXk5MDAHBycioztl69egBQYoh1WFiY0c9//PEHLl26hKCgIKOhxB4eHpg4caJRrFarxbfffou2bdti2LBhRucGDBiANm3a4MCBA4b5/6buWZo//vgD586dQ48ePQxD//UmT56MJk2aYP/+/aVuSVgepqYT7N+/H0IIDB482Oj4hAkT8P7775dYx8HPzw8ASh2aPWHCBKM52vqRHSqVCkDV9/W4ceMMw+Pt7OwMT6ojIyMhlUoBABKJBO3btwcA3L17t9KfS1l9Yan+MCUoKAj29vZG0wl0Oh2OHj2KPn36mNyZojJ5UZry/JuoU6cOlixZAo1GgyVLlmDBggWwt7fHihUrDJ9XWZRKJQDg+vXrRsevXbuG3bt3o1WrVhg7dmy5rkVERCVxKgEREVlc3bp1kZ6ejsLCwjIXCdOvpv700GpPT0+jn/XDpX19fUtco0OHDkY/JycnIy8vDxqNptT95XU6HQoKCpCcnGz4QlLaPUvz+++/A4DR8Hw9mUwGX19fHDp0CMnJyYbV1iuiZcuW8PX1RUxMDPLy8uDo6AiguDDQoEEDBAUFGcXrv8impqYiPj4et2/fRmJiIn755RcAxXPRn9a8eXOjn+vWrQugePs7oOr7ulmzZkY/6/vg6dfot31Uq9WV/lzK6gvAMv1hiqOjI3r27IkffvjBkA/nz5/HgwcP0K9fP5Ovq0xelKY8/yYAoFevXhg2bBh2794NoHiBwKf71hxTOxMsXboUOp0OCxYs4IKDRER/An+DEhGRxTVt2hTp6em4efNmmSum6/dO9/DwMDqu/9Knl52dDaD0UQj6UQd6WVlZAIDExESsXbvW5L311zR1z9LoR0PUqVOn1PP6eel/Zvu4wYMHY+nSpYiJiUH//v0RFxeH69evY+zYsZDJZEaxd+7cwZIlSxATEwMhBKRSKVq0aAFfX98ST1/1nt4KUf90XggBoOr72sHBoVztflJlP5ey+gKwTH+YExYWhmPHjuHUqVPo168fDh06BFdX11KLIHqVyYvSlOffxJPt1BcGOnXqVO7XAYCbmxvc3NyQkJAAIQQkEgn27duH8+fPo3///iXWHSAioophYYCIiCyud+/euHTpEo4fP262MFBUVISTJ0/CxsYGwcHBZq+pXxCttC9UT3/Z03+Be/HFF7F8+fIKtt48/bWfXg1fT//FsKztCc0ZMGAAVqxYgYMHD6J///7Yv38/AJTYjUAIgcmTJ+PGjRt45ZVXEBoaCoVCATs7OyQlJWHPnj2Vun916euKeJ6fS1X3R0hICGQyGY4cOYK+ffviyJEj6NOnj8lh+s8rL8wpKCjAsmXL4OjoCI1Gg3nz5iEqKqpCT/mVSiV+/PFHpKSkwNXVFStXroSjo2OZO4QQEVHZuMYAERFZXHh4OJycnLBt2zaTX9QAYPv27bh37x6CgoLQqFEjs9ds27YtJBIJLl26VOLclStXjH729vaGXC7H1atXS73Wp59+io8//hi5ubllv5mn6Ic8l7alnRACFy5cgKOjI5o0aVLha+u5uLggICAAp0+fRkFBAQ4dOgSFQlGiyBIXF4eEhAT069cPb775Jv7yl78YnvDqR2I8+eS7vKpLX1fE8/xcqro/nJ2d0a1bN5w6dQrnz5/HvXv3zE4jqGhemNrusCI+/PBD3Lp1C//+97/x0ksv4dq1a/j0008rdA391IqEhAR88sknSE1NxWuvvVbm7wYiIiobCwNERGRxjRo1wpw5c5CZmYkJEyYgOTm5RExUVBTee+89ODs7Y9GiRWVes2HDhujRowd++OEHnDhxwnD8wYMH2LRpk1GsnZ0d+vfvj/j4eGzZssXo3MmTJ/H+++/jwIED5Voc8WlNmjSBn58fLl68iK+//tro3McffwyVSoWwsLAKDcEuzeDBg5Gfn48tW7bg9u3bJRYdBP5/mPeDBw+MjqelpWHVqlUAjOfKl1d16euKeJ6fS3Xoj7CwMGRlZeG9995DgwYNDIsIlqaieaEfeaDVaivUJr2LFy9i69ataNeuHcaOHYtXX30VXl5eWLduHRITE8t9HX1h4NixY9i8eTO8vb0xfvz4SrWJiIiMcSoBERFViYiICGi1WixbtgwDBw5Er1690KpVKxQUFODs2bNISEiAl5cX1qxZgxdeeKFc11ywYAFGjhyJKVOmICwsDC4uLjh8+HCpc9Jnz56NCxcu4D//+Q+OHj0KX19f3L17F0ePHoWDgwOWLFlS6fe2ePFijBo1CnPnzsXBgwfRsmVLXL58GefPn0fz5s0xe/bsSl9b729/+xvq1auH9evXQyqVYuDAgSVivL290b59e5w9exZjxoxBx44d8fDhQxw9ehRA8aJ7GRkZlbp/denrinien0tV90doaCgWLVqEX3/9FREREWZX+69oXuj//e3YsQMPHjxAZGRkudulVqsxb948AMA777wDqVQKqVSKBQsWYNKkSZg3bx527txpcveEJ+lHfURHRwMo7vOn19QgIqLK4YgBIiKqMqNGjcKBAwcwYsQIqFQq7NixA/v374eTkxMWLlyIvXv3om3btuW+XosWLbBr1y706dMHZ86cwbfffouePXuWOpfbxcUFX331FcaPH4+7d+9i27ZtuHTpEvr06YOvvvoKHTt2rPT78vb2RnR0NF588UX8/vvv2L59O9LS0jBp0iRER0fDxcWl0tfWk8vl+Pvf/478/Hz4+/sbFs97kkQiwbp16zB48GDcunUL27Ztw4ULFxAaGoo9e/agY8eOuHLlSqWKA9WlryvieX4uVd0frq6u6Ny5MwCgb9++ZmMrmhddu3bF6NGjkZ6ejh07dlToKf/q1auRlJSEUaNGGbaPBIq3Wezbty9+/fVXfP755+W6VsuWLQ2FgL59+6Jnz57lbgcREZknEZWZXEhEREREREREtQJHDBARERERERFZMRYGiIiIiIiIiKwYCwNEREREREREVoyFASIiIiIiIiIrxsIAERERERERkRVjYYCIiIiIiIjIirEwQERERERERGTFWBggIiIiIiIismIsDBARERERERFZMRYGiIiIiIiIiKwYCwNEREREREREVoyFASIiIiIiIiIr9n8dja+9CLKuLAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=[12, 5])\n", "ax.set_yscale('log')\n", "data_smooth.plot(ax=ax)\n", "ax.set(\n", " title='Ill-Conditioning of Vandermonde Matrix $V$', \n", " xlabel='Order of Vandermonde Matrix $V$', \n", " ylabel='Value', \n", " yscale='log', \n", " yticks=[1,1e100,1e200]);" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.7 ('base')", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" }, "vscode": { "interpreter": { "hash": "ad2bdc8ecc057115af97d19610ffacc2b4e99fae6737bb82f5d7fb13d2f2c186" } } }, "nbformat": 4, "nbformat_minor": 1 }